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1. Introduction 

Final states containing hadronic jets are produced at large rates at high energy particle 
colliders. Since the distribution of these jets relates directly to the dynamics of the hard 
parton-level scattering process, jet observables can be used to probe elementary particle 
reactions. Owing to their large production cross sections, various jet observables can be 
measured to a high statistical accuracy, thus making them ideal for precision studies [1]: 
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for example three-jet production in e^e~ annihilation to determine the strong coupling 
constant ag and for quantum chromodynamics (QCD) studies, two-plus-one-jet produc- 
tion in deep inelastic scattering as a probe of the gluon distribution in the proton and to 
measure a^, single jet inclusive production and vector-boson-plus-jet production at hadron 
colliders as benchmark reactions, and to measure parton distributions. Experimental data 
on these observables often attain an accuracy of a few per cent or better, such that mean- 
ingful precision studies must rely on theoretical predictions accurate to the same level. 
In perturbative QCD, this precision usually requires corrections at next-to-next-to-leading 
order (NNLO). 

NNLO calculations of observables with n jets in the final state require several ingre- 
dients: the two-loop corrected n-parton matrix elements, the one-loop corrected (n -|- 1)- 

parton matrix elements, and the tree-level (n -|- 2)-parton matrix elements. For most 
massless jet observables of phenomenological interest, the two-loop matrix elements have 
been computed some time ago [2, 3], while the other two types of matrix elements are 
usually known from calculations of next-to-leading order (NLO) corrections to {n + 1) jet 
production [4,5]. 

The one-loop {n + l)-parton matrix elements contribute to n jet observables at NNLO 
if one of the partons involved becomes unresolved (soft or collinear) [6]. In these cases, 
the infrared singular parts of the matrix elements need to be extracted and integrated 
over the phase space appropriate to the unresolved configuration to make the infrared pole 
structure explicit. Methods for the extraction of soft and collinear limits of one-loop matrix 
elements are worked out in detail in the literature [6-10]. Likewise, the {n + 2)-parton 
matrix elements contain double real radiation singularities corresponding to two partons 
becoming simultaneously soft and/or collinear [11-14]. To determine the contribution to 
NNLO jet observables from these configurations, one has to find two-parton subtraction 
terms which coincide with the full matrix element and are still sufficiently simple to be 
integrated analytically in order to cancel their infrared pole structure with the two-loop 
virtual and the one-loop single-unresolved contributions. Often starting from systematic 
methods for subtraction at NLO [15-18], several NNLO subtraction methods have been 
proposed in the literature [19-23], and are worked out to a varying level of sophistication. 

For observables with partons only in the final state, an NNLO subtraction formalism, 
antenna subtraction, has been derived in [24]. The antenna subtraction formalism con- 
structs the subtraction terms from antenna functions. Each antenna function encapsulates 
all singular limits due to the emission of one or two unresolved partons between two colour- 
connected hard radiator partons. This construction exploits the universal factorization of 
matrix elements and phase space in all unresolved limits. The antenna functions are de- 
rived systematically from physical matrix elements [25]. This formalism has been applied 
in the derivation of NNLO corrections to three-jet production in electron-positron annihi- 
lation [26,27] and related event shapes [28,29], which were used subsequently in precision 
determinations of the strong coupling constant [30-34]. The formalism can be extended 
to include parton showers at higher orders [35], thereby ofi'ering a process-independent 
matching of fixed-order calculations and logarithmic resummations [31,36], which is done 
on a case-by-case basis for individual observables [37] up to now. The formalism can be 
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extended to include massive fermions [38]. For processes with initial-state partons, antenna 
subtraction has been fully worked out only to NLO so far [39]. In this case, one encoun- 
ters two new types of antenna functions, initial-final antenna function with one radiator 
parton in the initial state, and initial-initial antenna functions with both radiator partons 
in the initial state. In this work, we derive all NNLO initial-final antenna functions and 
perform their integration over the appropriately factorized phase space. These functions 
form part of the full set of antenna functions needed for NNLO calculations of hadron 
collider processes, and are, together with the already known final-final antenna subtraction 
terms, sufficient for NNLO calculations of jet observables in deeply inelastic lepton-hadron 
scattering. 

Other approaches to perform NNLO calculations of exclusive observables with initial 

state partons are the use of sector decomposition and a subtraction method based on the 
transverse momentum structure of the final state. The sector decomposition algorithm [40] 
analytically decomposes both phase space and loop integrals into their Laurent expan- 
sion in dimensional regularization, and performs a subsequent numerical computation of 
the coefficients of this expansion. Using this formalism, NNLO results were obtained for 
Higgs production [41] and vector boson prodiiction [42] at hadron colliders. Both reactions 
were equally computed independently using an NNLO subtraction formalism exploiting 
the specific transverse momentum structure of these observables [43]. 

This paper is structured as follows: in Section 2, we construct the subtraction terms 
required at NNLO for initial-final configurations with two unresolved partons at tree-level or 
one unresolved parton at one loop. They require two new types of antenna functions: tree- 
level 2^3 and one-loop 2 — > 2, each with one parton and one off-shell neutral current in 
the initial state. The phase space mappings relevant to these antenna subtraction terms are 
discussed in Section 3. The analytic integration of the new initial-final antenna functions is 
described in Section 4, it proceeds through a reduction to a set of master integrals, which 
we computed and collect in Appendix A. The integrated antenna functions are listed in 
Sections 5 and 6 for the quark-initiated and gluon-initiated antennae. A strong check 
on these results is provided by the rederivation of NNLO coefficient functions for deep 
inelastic scattering, which we carry out in Section 7. Finally, we conclude with an outlook 
in Section 8. 



2. Initial-final antenna subtraction at NNLO 



Antenna subtraction of initial-final configurations at NLO is derived in detail in [39] . Sub- 
traction terms in the case of one hard parton in the initial state are built in the same 
fashion as for the final-final case (formula (2.5) in [24]). Wc have the following subtraction 
term associated to a hard radiator parton i with momentum p in the initial state: 

da^'^'-f\p, r) = V d^rn+iih, • • • , km+i;p, r) 

m+l 

X ^ Xo^.fe \Mm{ki, ...,Kk,..., km+i;xp, r)|2 jM {ki,...,KK,..., km+l) . (2.1) 
j 
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The additional momentum r stands for the momentum of the second incoming particle, 
for example, a virtual boson in DIS, or a second incoming parton in a hadronic collision 
process. This contribution has to be appropriately convoluted with the parton distribution 
function /j. The tree antenna -^^ depending only on the original momenta p, kj and 
kk, contains all the configurations in which parton j becomes unresolved. The m-parton 
amplitude depends only on redefined on-shell momenta fci , . . . , Kk , ■ ■ ■ , and on the mo- 
mentum fraction x. In the case where the second incoming particle is a parton, there 
is an additional convolution with the parton distribution of parton r and corresponding 
subtraction terms associated with it. 

The jet function, in (2.1) depends on the momenta kj and kk only through Kk- 

Thus, provided a suitable factorization of the phase space, one can perform the integration 
of the antennae analytically. Due to the hard particle in the initial state, the factorization 
of phase space is not as straightforward as for final-final antennae. 

The phase space can be factorized in an m-parton phase space convoluted with a two 
particle phase space: 



d^m+i (ki,..., km+i ■,P,r) = d^rn (ki, . . . , Kk , ■ ■ ■ , km+i ; xp, r) 

—d^2{kj,kk;p, q) — 

ZTT X 



X ^d^2{kj,kk;p,q)—, (2.2) 



where = —q^ = —{kj 4- fcfc ~ p)'^- Replacing the phase space in (2.1), we can explicitly 
carry out the integration of the antenna factors over the two particle phase space. When 
combining the integrated subtraction terms with virtual contributions and mass factoriza- 
tion terms, it turns out to be convenient to normalize the integrated antennae as follows 

where 

C{e) = {A^r . (2.4) 

The integrated form of the subtraction term is then 

da^'^'f\p,r) = EE/^ / — X^,jkM^)d^m{h,...,KK,...,km+i;xp,r) 

m+l j ^ 

X \Mm{ki, Kk, km+i;xp, r)f jjn\ki, Kk, km+i){2.5) 

Finally, the subtraction term has to be convoluted with the parton distribution functions 
to give the corresponding contribution to the hadronic cross section. The explicit poles 
in the integrated form cancel the corresponding ones in the virtual and mass factorization 
contributions. To carry out the explicit cancellation of poles, it is convenient to recast, 
by a simple change of variables, the integrated subtraction term, once convoluted with the 
parton distribution functions (PDFs), in the following form 
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(2.6) 



This convolution has already the appropriate structure and mass factorization can be car- 
ried out explicitly leaving a finite contribution. The remaining phase space integration, 
implicit in the Born cross section, do"^, and the convolutions can be safely done numeri- 
cally. When considering reactions with only one incoming hadron, the second PDF has to 
be replaced by a Dirac delta. Reactions with two hadrons will require additional subtrac- 
tions containing initial-final antennae involving the second parton in the initial state and 
initial- initial antennae as well. This case is discussed to NLO in [39]. 

At NNLO, two types of contributions to m-jet observables require subtraction: the tree- 
level m + 2 parton matrix elements (where one or two partons can become unresolved), and 
the one-loop m + 1 parton matrix elements (where one parton can become unresolved). 
The corresponding subtraction terms are denoted by daffp^^Q and da^f^^Q. Final-final 
antenna subtraction terms for both cases are constructed in [24]. 

In dcr^jY^Q, we have to distinguish four different types of unresolved configurations: 

(a) One unresolved parton but the experimental observable selects only m jets; 

(b) Two colour-connected unresolved partons (colour-connected); 

(c) Two unresolved partons that are not colour connected but share a common radiator 
(almost colour-unconnected); 

(d) Two unresolved partons that are well separated from each other in the colour chain 
(colour- unconnected) . 

Among those, configuration (a) is properly accounted for by a single tree-level three-parton 
antenna function like used already at NLO. Configuration (b) requires a tree- level four- 
parton antenna function (two unresolved partons emitted between a pair of hard partons) , 
while (c) and (d) are accounted for by products of two tree-level three-parton antenna 
functions. With radiator parton i in the initial state, the subtraction terms for these 
configurations read: 




(2.7) 
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X \Mm{ki, Kl, km+2;xp, r)^ jjrP^ki, Kl, km+2) 



(2.8) 



^'^NNLO = y] d$m+2 {kl, ■ ■ ■ , km+2 P, r) — ^ 

m+2 ^"^+2 

X] -^hk XmiK \Mm{ki, Kk, Km, km+2;xp, r)|" 

. j,l 

Xj^\ki,...,KK,KM,...,km+2) 

+ X] ^kim ^IjK \Mm{ki, Kk, Km, km+2);xp, r\ 

hi 

xJ^\ki, ...,Kk, Km, km+2) 

^(^NNLO = ^<^^m+2{ki,. . . ,km+2;V,r)-^ 



m+2 



'm+2 



X 



X -^hk x^i,im \Mm{ki, Kk, Km, km+2;xp, r)|^ 



3,1 



xJ^\ki,..., Kk, Km, • • • , km+2) 
+ X ^%jk \Mm{ki, Kk, Km, km+2);xp, \ 

0,1 

xJ^\ki,...,KK,KM,...,km+2) 
^(^NNLO = ^<^^m+2{ki,...,km+2;P,r)-^ 



m+2 



'm+2 



X 



(2.9) 



(2.10) 



X -^hk ^nap \Mm{k\, • • • , Kk, • • • , Kn, Kp, km+2;xp, r)f 



3,0 



X J^\ki, . . .,Kk, . . .,Kn,Kp, km+2) 



(2.11) 



As before, the original momenta of the (m + 2)-parton phase space are denoted by j,k, . . ., 
while the combined momenta obtained from a phase space mapping are labelled hy J, K, . . .. 
Only the combined momenta appear in the jet function, x^j^^ is a three-parton sub-antenna 
function containing only limits where parton b is unresolved with respect to parton a, but 
not limits where parton b is unresolved with respect to parton c. dcr^^^^ applies if the 
common radiator is in the final state, while da^^^^ applies if the common radiator is in the 
initial state. The only genuinely new ingredient here is the four-parton initial-final antenna 
function X^ji^i, which can be obtained by crossing the corresponding final- final antenna 
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functions, and has to be integrated analytically over the appropriate antenna phase space. 
The resulting integrated antenna function exploits the phase space factorization 



d^m+2{ki, km+2;p, r) = d^rn{ki, ...,Kl,..., km+2;xp, r) 



X ^d<t2,{kj,kk,ki;p,q)— , 

ZTT X 



yielding the integrated antenna 



■yO _ 



'27r' 



(2.12) 



(2.13) 



In all products of two three-parton antenna functions, the analytic integration has to 
be performed only over the outmost antenna function, yielding the integrated NLO antenna 
functions in the appropriate (final-final or initial-final) kinematics. 

The one-loop single unresolved subtraction term da^^]^Q must account for three types 
of singular contributions: 



(a) Explicit infrared poles of the virtual one-loop (m -|- 1) parton matrix element. 

(b) Single unresolved limits of the virtual one-loop (m -|- 1) parton matrix element. 

(c) Terms common to both above contributions, which are oversubtr acted. 

With radiator parton i in the initial state, the subtraction terms for these configurations 
read: 



AVS,l,a 

'^"nnlo 



J\f d^m+i{ki, . . . ,km+i;p,r)—^ 

^ ~ J ~x . . . , fcjfc, . . . , km+i;xp,r) 

ik 

XJ^+'\ki,...,kk,...,km+l) 

^(^NNLO = d$m+i(fci, • • ■ ,km+i;p,r)—^ 

Drr ' 



(2.14) 



m+l 



J2 Kjk \Ml,{ki, ...,Kk,..., km+i;xp, r)|2 J^\ki, ...,Kk,..., km+i) 

3 - 

+ ^Ijk \Mm{ki, . . . , Kk, km+i;xp, r)p J^^ {ki, Kk, km+i) 



(2.15) 



J ~x ^ ^""P \Mmiki, ...,kk,..., Kn, Kp, km+i;xp, r)p 
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xJ^\ki, ...,kk,..., Kn, Kp,..., km+i) 



(2.16) 



m+l 



m+l 



X '^nopi^np) X l-^m(A;i, ■ ■ ■ , Kk, . . . , kn, kp, . . . , km+i;xp, r)\' 



np 



xjj^\ki, Kk, ...,kn,kk,..., km+i) 



(2.17) 



In here, ^Ij^ denotes a one-loop three-parton initial-final antenna function, which is the 
only new ingredient. These antenna functions can be obtained by crossing from their final- 
final counterparts, listed in [24], and have to be integrated over the appropriate phase 
space: 

J^/^^gx/,.. (2.18) 



The subtraction terms dcr^^^Q, do'ffNLO ^^'^ ^'^^m.o require three different types of 
antenna functions corresponding to the different pairs of hard partons forming the antenna: 
quark-antiquark, quark-gluon and gluon-gluon antenna functions. We derived these an- 
tenna functions [25] for final-final kinematics in a systematic manner from physical matrix 
elements known to possess the correct limits. For the initial-final kinematics, one parton is 
crossed into the initial state. Special care has to be taken in these crossings, if we start from 
those final-final antenna functions which contain more than one quark/ antiquark or more 
than one gluon in the final state. In the case of more than one quark/antiquark pair of 
different fiavour (final-final antenna functions B^, E^, £'3 -type, £^4 -type), we have to distin- 
guish the crossing of the primary quark q (which is coupled to the external current) and the 
secondary quark q' (which is not coupled to the external current). The identical flavour 
antenna function C2{lq,2q,Aq,2q) is constructed from the interference of the four-quark 
amplitudes with the antiquark-momenta interchanged, and it contains only the (3^, 4q, 2q) 
triple collinear limit. Consequently, it is symmetric in the two antiquark momenta, but not 
in the two quark momenta, and has thus three different crossings: Cqqqq (either antiquark 
crossed), qqq (quark (1^) crossed) and C^qqq (quark (3^) crossed). Crossing one of several 
gluons into the initial state is unambiguous for most antenna functions owing to their sym- 
metry properties. The only exception is the quark-gluon antenna function .04(15, 3g,Ag,5g), 
where gluons (3g) and (5^) are colour-connected to the quark, while gluon (4^) is not. We 
thus distinguish two crossings, D^^^^ (gluon (3^) or (5g) crossed) and I?^, (gluon (4^) 



crossed). We list all NLO and NNLO initial- final antenna functions with an initial state 
quark in Table 1 and with an initial state gluon in Table 2. 

It was shown in [27] that these antenna subtraction terms result in an oversubtraction 
of large-angle soft gluon radiation. To correct for this oversubtraction, one introduces the 
soft antenna function 

(2.19) 
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Quark initiated tree level 



one loop 



quark-quark 










A'^ 

9,99' 9,99' 9,99 


q gm 


40 




q ^ q'q'q 


^q,q'q'q 




q' qqq' 


^q',qqq' 




q qqq 


(^0 (^0 /^O 

q,qqq^ qmq^ qmq 





quark-gluon 






Q^gg 


q,99 


9,99' 9,99 


Q ggg 


q,999 




q q'q' 


^q,q'q' 


^q,q'q" ^q,q'q" ^q,q'q' 


q q'q'g 


q,q'q'9' q,q'q'9 




q' q'q 


^q',q'q 


^q',q'q' ^q',q'q' ^q',q'q 


q' q'qg 


q',q'q9' q',q'q9 




gluon-gluon 






Q^Qg 


q,q9 


f^i fii f^i 
q,q9' 9,99' 9,99 


Q qgg 


^q,q99' ^q,q99 




q qq'q' 


^q,qq'9' 





Table 1: List of tree level and one loop antenna functions for the initial-final configurations with 
a quark in the initial state. 



where a and c label arbitrary hard partons. Those soft factors are associated with an 
antenna phase space mapping {i,j,k) iI,K) (final-final) or {p,j,k) {xp,K) (initial- 
final). In contrast to all previous subtraction terms, the hard momenta a, c do not need 
to be equal to the hard momenta of partons i, k in the antenna phase space - they can be 
arbitrary on-shell momenta in the initial or final state. 

If parton (a) is in the initial state, and (c, i, j, k) are in the final state, the integral of 
each of the soft antenna function over the antenna phase space can be written as 



(sik) 



d^XijkSajc 

r2(l-e)e^^ / 2 



r(l-3e) V e 



where we have defined 



XacJK 



+ In {xac,iK) + eLi2 



tacSlK 



1 — X, 



ac,IK 



XacJK 



{tal +taK)iScI + Sck) 



(2.20) 
(2.21) 
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Gluon initiated tree level 



one loop 



quark-quark 



g^qq 
9 gqq 


gm 

M /I 

fl>fl99' gmg 


2i 1 ^1 

ff>9?' fl>99' gm 


quark-gluon 






9^ 99 
9 99q 

9 qq'q' 


gm 

^gmq' g',ggq 
gm'g''' gm'g' 


g,gq^ gm 


gluon-gluon 






9^ 99 


9,9g 


pi pi 
g,gg^ g,gg 


9 999 

9^qq 
9 qq9 


po 
g,ggg 

gm 

gmg^ gmg 


/^i f^i Ai 
gmi gm^ gm 



Table 2: List of tree level and one loop antenna functions for the initial-final configurations with 
a gluon in the initial state. 



If parton (i) is in the initial state, while (o, c, j, k) are in the final state, we obtain the 
following integral: 



where we have defined 



""'^^ (SaK + (1 - x)Sai) (ScK + (1 - x)Sci) ' 



3. Phase space mappings 

As discussed above, the construction of subtraction terms requires mapping the original 
set of momenta onto a reduced set. The mappings interpolate between the different soft 
and collinear limits which the subtraction term regulates. Appropriate mappings for the 
initial-final configurations, both for single and double unresolved configurations have been 
discussed in [39], and are only briefly summarized here. 

The proper subtraction of infrared singularities requires that the momentum mapping 
satisfy 

xp ^ p , Kk — > fefc when j becomes soft , 
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xp ^ p , Kk kj + kk when j becomes coUinear with k , 

xp ^ p — kj , Kk kk when j becomes coUinear with i . (3-1) 

In this way, infrared singularities are subtracted locally, except for angular correlations, 
before convoluting with the parton distributions. That is, matrix elements and subtraction 
terms are convoluted together with PDFs. In addition, the redefined momentum, Kk, 
must be on shell and momentum must be conserved, p — kj — kk = xp — Kk, for the phase 
space to factorize as above. This is accomplished by: 

_ Sij + Sik - Sjk 

X — , 

Slj + Sik 

Kk = kj + kk-{l-x)p, (3.2) 

where sij = {p — kj)'^, etc. If parton j becomes soft or coUinear to parton A;, x ^ 1. If 
parton j becomes collinear with the initial state parton i, x = 1 — z with z the fraction of 
the momentum p carried by parton j. 

This mapping is easily generalized to deal with more than one parton becoming un- 
resolved. As explained above, the building blocks for the double real radiation in the 
initial-final situation are colour-ordered four-parton antenna functions Xij^i, with one ra- 
diator parton i (with momentum p) in the initial state, two unresolved partons j, k and 
one radiator parton / in the final state: 

^ _ Slj + Slk + SU — Sjk — Sjl — Ski 

Slj + Slk + Sll 

Kl = kj + kk + ki- {1- x)p , (3.3) 

where kj, kk and ki are the three final-state momenta involved in the subtraction term. It 
satisfies the appropriate limits in all double singular configurations: 

1. j and k soft: x — > 1, K^ — > ki, 

2. j soft and kk \\ kf x — >■ 1, Kl kk + ki, 

3. kj = zp \\ p and kk soft: x — > 1 — z, Kl — > h, 

4. kj = zp \\ p and kk \\ kf. x ^ 1 — z, Kl kk + ki, 

5. kj II kk II kf. X 1, Kl — > kj + kk + h, 

6. kj + kk = zp \\ p: X ^ 1 — z, Kl — >■ h, 

where partons j and k can be interchanged in all cases. 

The construction of NNLO antenna subtraction terms requires moreover that all single 
unresolved limits of the four-parton antenna function X^jkl have to be subtracted, (2.8), 
such that the resulting subtraction term is active only in its double unresolved limits. A 
systematic subtraction of these single unresolved limits by products of two three-parton 
antenna functions can be performed only if the NNLO phase space mapping turns into an 
NLO phase space mapping in its single unresolved limits. 

In the limits where parton j becomes unresolved, wc denote the parameters of the 
reduced NLO phase space mapping (3.2) by x' and K'j^. We find for (3.3): 
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1. j becomes soft: 



X — — = x' , Kl ^ kk + k - {1 - x)p = K'j^ . 



2. kj II kk, kj + kk = Kk- 



X 



SlK + Sil - SKI 
SlK + Sll 



= X 



KL^kK + ki-{l- x)p = K'l 



3. kj = zp II p: 



sik + su 



KL^kk + ki-{l- x'){l - z)p = K'l 



It can be seen that in the first two hmits, the NLO mapping involves the original incoming 
momentum p, while in the last limit (initial state collinear emission), it involves the rescaled 
incoming momentum (1 — z)p. To subtract all three single unresolved limits of parton j 
between emitter partons i and k from X^jkh one needs to subtract from it the product 
of two three-parton antenna functions Xijk • Xi^ki- The phase space mapping relevant 
to these terms is the iteration of two NLO phase space mappings. Analytical integration 
of terms with this mapping is required only over the phase space appropriate to the first 
antenna function. 

Equally, parton k can become unresolved. Expressing the reduced NLO phase space 
mapping by x" and K'^. We find for (3.3): 



1. k becomes soft: 



'Ji±I^^ = :,"^ KL^kj + k-{l-x)p = Kl 
sij + su 



2. kk II kj, kj + kk = Kk- 



SlK + Sll - Ski 

SlK + Sll 



x". 



KL^kK + ki-(l- x)p = K'l . 



3. kk II k, ki + kk = Kk: 

SlK + Sij - SKj 



X 



SlK + Sij 



X 



KL^kK + kj - (1 - x)p = K'l . 



In all limits, the reduced NLO mapping involves the original incoming momentum p. Con- 
sequently, the three single unresolved limits of parton k between emitter partons j and I 
can be subtracted from Xijki by a product of a final-final and an initial-final three-parton 
antenna function Xj^r Xi^j^. The phase space mapping relevant to these terms is the prod- 
uct of an NLO final-final phase space mapping with an initial-final mapping. Integration 
of the final-final antenna phase space yields a constant, not involving an extra convolution. 
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4. Integration of initial-final antenna functions at NNLO 



The initial-final antenna functions all have the scattering kinematics 

q + Pi-'Pi +P2{+P3) , 

where 

2 /n2 ^ r\ 2n ^ 222n 
q =-Q <0, p- =0,z = , p^=p^=p^ = 0, 

^Q-Pi 

and p3 is present only for the NNLO real radiation antenna functions. Consequently, 
integration over the final-state two-parton or three-parton phase space yields a result which 
depends only on and z. From dimensional counting, one can immediately conclude 
that the dependence on is only multiplicative, according to the mass dimension of the 
integral. 

The NNLO double real radiation antenna functions X^jf^^ have to be integrated over 
the inclusive three-parton final state phase space. The NNLO one-loop single real radiation 
antenna functions X^j}. are integrated over the inclusive two-parton final state phase space, 
and over the loop momentum. For both types of integration, wc employ the by-now stan- 
dard technique of reduction to master integrals. The master integrals are then computed 
from their differential equations. 

To perform the reduction, we first express all phase space integrals as loop integrals 
with cut propagators [44] . Consequently, all NNLO integrals are expressed as cuts of two- 
loop four-point functions with two off-shell legs in forward scattering kinematics. Using 
integration-by-parts (IBP, [45]) and Lorentz invariance (LI, [46]) identities among the in- 
tegrals of any given topology, the large number of different integrals can be expressed in 
terms of a small number of master- integrals. This reduction is performed iteratively, based 
on the lexicographic ordering of the integrals, expressed by the Laporta algorithm [47]. 

After carrying out the reduction, one finds nine master integrals for the NNLO double 
real radiation antenna functions, described in Section 4.1 below, and listed in Appendix A.l. 
For the NNLO one-loop single real radiation antennae, one finds six master integrals, which 
we describe in Section 4.2 and list in Appendix A. 2. 

We computed these integrals both directly and by using the differential equation tech- 
nique [46,48,49]. To derive the differential equations for each master integral, we employ 

to carry out the differentiations at the integrand level. The boundary conditions required 
for the solution of the differential equations are either obtained from self-consistency con- 
ditions on the integrals, or by explicit evaluation. The explicit evaluation is very similar 
to the evaluation of inclusive four-point phase space integrals, described in [50]. 

Some of the more involved master integrals can be related to phase space integrals 
computed by Zijlstra and van Neerven in the context of the NNLO corrections to deep 
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inelastic structure functions [51]. Where appropriate, we compared our results to the 
expressions in the appendix of [51], finding full agreement. Explicit expressions for all 
master integrals are listed in Appendix A.1-A.2. 

All master integrals contain multiplicative factors of the form (1 — z)"^ or (1 — z)~'^^, 
which regulate soft endpoint singularities in initial state convolution integrals. These fac- 
tors must be left as such in the master integrals, and can be expanded in the form of 
distributions 

(1 _ = _ ^(1 - + ^ -ir-^n{z) , (4.3) 



n 

with 



ln"(l 



l-z 

only after being inserted into the integrated antenna functions. All other terms in the mas- 
ter integrals can be expanded, yielding Harmonic Polylogarithms (HPLs, [52]) of argument 
z. 

4.1 Tree-level 2^3 antenna functions 

For the tree-level 2 — > 3 antenna functions, we have the 'DIS-like' process 

q + Pi^Pi+P2+P3 ■ 

There arc 12 propagators, including the three that are cut in the phase space integration 
(L>io, Dn, Du): 





= {Pi-Pif 


D2 


= {Q-Pif, 




= {P2 +P3)^ 


Da 


= {Pi - P2f 




= {Q- P2f , 


De 


= {Pl +P3f 


Dr 


= {Pi - Psf 


Ds 


= {q- p^f , 


Dg 


= {Pl +P2f 


Dio 


= pI, 


Du 


= pI , 


Dl2 


=ph 



(4.4) 

To perform the reduction to master integrals, we impose momentum conservation = 
q + Pi — Pl — P2, set p? = 0, q^ = —Q^ and drop any integral where Dio, Du and D12 
are not in the denominator. After labelling the inclusive phase space integral as I[0], the 
convention for naming the master integrals follows the labelling of the numerators, i.e. 

/[1. 2, 51 = / MJMjM + „ - - - , (4.5) 

J J^2 J^b 
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I[0] I[2] I[2,6] 




I[l,2,5] I[2,3,5] I[2,4,9] 




I[l,3,4,6] I[2,3,5,6] I[l,2,4,5] 



Figure 1: Master integrals for the phase space integration of the tree level initial- final antennae at 
NNLO. The double line in the external states represents the off-shell momentum, q with = —Q^, 
the single one is the incoming parton. All internal lines are massless. The cut propagators are the 
ones intersected by the dotted line. 



master 


deepest pole 


behaviour at a: = 1 


known to order 


m 






all 


I[2] 




(1 - 


all, 


I[2,6] 




(1 - 


all, 


7[1,2,5] 


e-2 


(1 - 




7[2,3,5] 


e-2 




e3 


/[2,4,9] 


e-3 


(1 - 


e3 


I[l,3,4,6] 


e-3 


(1 - .t)-1-2^ 


all, 


I[2,3,5,6] 


€-3 


(1 -a;)-i-2. 


el 


I[l,2,4,5] 


e-2 


(1 - 





Table 3: Summary of the main properties of the three particles phase space master integrals 
where 

When squaring the 2^3 antennae, we find at most 4 propagators, plus the 3 cut ones. 
All the integrals can be reduced to the set of 9 master integrals shown in Figure 1. All the 
masters, except /[1, 2, 4, 5], have been computed by direct integration and by the differential 
equations method, sTipplemented, where necessary, by a direct calculation at a; = 1 after 
factorizing the leading singularity. The x = I boundary conditions for /[1,3, 4, 6] and 
/[2, 3, 5, 6] were checked numerically using sector decomposition [40]. /[I, 2, 4, 5] has been 
computed only using the differential equations method. The master integrals /[I, 2, 5], 
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Topology 1 


Topology 2 


Topology 3 


Topology 4 




K 


K 


1.2 


1.2 


Do 


(h 4- r\ \'^ 


\K -t- Pi) 




(h -L- n \'^ 
[K -\- Pi) 


Ds 




{k+pi-pif 


{k + pi-pxf 


{k+pi-pif 


Di 


{k-q? 


{k-qf 


{k+P2f 


{k+P2f 


D5 


ipi - Pif 


{q-pi? 


{Pi-Plf 


{q-pi? 


Da 


p\ 


pI 


pI 


pI 


D7 


pI 


pI 


pI 


pI 



Table 4: Definition of the topologies for the combined phase space and loop integration. 




V[l,3] V[l,41 V[2,4] 




V[l,3,4] V[l,2,3,4,5] C[l,2,3,4] 



Figure 2: Master integrals for the loop plus phase space integration of the one loop initial-final 
antennae at NNLO. The double line in the external states represents the off-shell momentum, q 
with = — Q^, the single one is the incoming parton. All internal lines are massless. The cut 
propagators are the ones intersected by the dotted line. 

/[2,3,5], J[2,3,5,6] and J[l,2,4,5] agree up to order e° with the results in [51]. We 
summarize in Table 3 some of the properties of the master integrals. 

4.2 One-loop 2^2 antenna functions 

When interfering the one loop 2^2 antennae with the corresponding tree level ones, we 
can combine the loop and phase space integrations. The partonic process in this case is 

q + Pi^Pi+P2 ■ 

Denoting the loop momentum by A;, we can identify four topologies, two planar (Topology 
1 and 2) and two non-planar (Topology 3 and 4). The topologies arc defined in Table 4. 
Topologies 1 and 2 only differ in the propagator Ds, as is also the case for Topologies 3 and 
4. Subtopologies of the non-planar integrals can be expressed by the planar topologies. 

All the resulting integrals can again be reduced to a small set of masters. In this case 
we find only six of them. They are shown in Figure 2. The notation for the master integrals 
follows the propagator definitions of Topology 2 in the planar case (denoted by V) and of 
Topology 3 in the non-planar case (denoted by C). 

The one loop plus phase-space master integrals have been computed using differential 
equations in external invariants, together with a direct calculation at a; = 1 after factorizing 



- 16 - 



master 


deepest pole 


behaviour at a; = 1 


known to order 


V [I, 6\ 


t 


[L X) 


all 

all, t 


V[1,A] 


e-i 


(I-X)- 


all, 


V[2A] 


e-i 


(1 - x)-^' 


all, 


F[l,3,4] 


e-i 


(l-x)-' 


e3 


^[1,2,3,4,5] 


e-3 


(1 -x)-l-2^ 


e3 


C[l,2,3,4] 


e-3 


(1 -x)-^ 


e3 



Table 5: Summary of the main properties of the one loop plus two particles phase space master 
integrals 



the leading singularity at x = 1. Integrals ^[1,2,3,4], y[l,2,3,4, 5] and C[l,2,3, 4] have 
been checked, up to order e*^, against a direct analytic calculation. We summarize in table 5 
some of the properties of the master integrals. Explicit expressions for the master integrals 
are listed in Appendix A. 2. 

We always integrate the unrenormalized 2—^2 one-loop squared matrix elements, 
divided by a normalization factor C(e), relevant to a particular antenna function, which we 
denote as X^'^^. The antenna function is obtained after renormalization and subtraction of 
the corresponding trcc-lcvcl antenna function multiplied with the one-loop correction to the 
hard radiator pair. Renormalization of the one-loop antenna functions is always carried out 
in the MS-scheme at fixed renormalization scale /x^ = Q^. It amounts to a renormalization 
of the strong coupling constant and (in the case of the quark-gluon and gluon-gluon antenna 
functions) to a renormalization of the effective operators used to couple an external current 
to the partonic radiators. The relation between renormalized and unrenormalized one- loop 
squared matrix elements is as follows: 

^l,R _ yl,U bo % yO (A 

^ = ^k^ (4-8) 

^- ;? = 4g - ^ ^ik - ^ , (4.9) 

where 

= ^ , 6o,F = A; (4.10) 
6 3 

are the colour-ordered coefficients of the one-loop QCD /3-function: 

/3o = 6oiV + ho,F^F . (4.11) 
The renormalization constants for the effective operators are 



?7o = , ?7o,F = for X = ^ , 
3 

/o = feo + 2 ' ^o,F = bo,F for X = D,E 

?7o = 2 6o , r/o,F = 2 bo p for X = F,G . 
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The one-loop antenna functions are obtained from the renormahzcd one-loop squared 
matrix elements by subtracting from them the product of the tree-level antenna function 
with the the virtual one-loop hard radiator vertex correction [6,24]: 





- ^i,jk 


-Rc(( 


-1)- 




i,jk 1 


(4.12) 




yl,R 
- ^i,jk 


-Re(( 


-!)■ 




Y^ 
i,jk ' 


(4.13) 




- Y^'^ 

- ^i,jk 


-Re(( 


-!)■ 




Y^ 


(4.14) 



The one-loop corrections to the hard radiator vertex are listed in [24] for Al, V^, V^, T\ 
and T\. Prom these, the remaining functions follow: 







^2 = 0, 






(4.15) 


v\ 


= 0, 








(4.16) 


El 


= vl 


4' = 0, 


si 


— ^2 ! 


(4.17) 




= 0, 








(4.18) 


Ql 


— -^2 


Ql = ^. 


Ql 


— -^2 • 


(4.19) 



5. Qucirk initiated antennae 

The quark-initiated initial-final antenna functions are obtained from the final-final antenna 
functions by crossing a quark into the initial state. Their unintegrated analytical expres- 
sions are obtained by a pure kinematical crossing from the final-final expressions listed 
in [24], with no extra symmetry factors or overall signs applied, and no further decompo- 
sitions into sub-antennae. This procedure differs from the approach used in our previous 
work on initial-final antenna functions at NLO [39]. In the application of the antenna 
subtraction formalism to the calculation of NNLO corrections to e~^e~ — > 3j [26], it turned 
out that several decompositions can be circumvented by an appropriate symmetrization 
of the real radiation matrix elements in the process under consideration, and that the re- 
maining required decompositions can be derived on a case-by-case basis, often not along 
the lines of [24,39]. A systematic decomposition would start from a specific process under 
consideration, and is beyond the scope of the present study. 

The virtual one-loop antenna functions are obtained from the final-final one-loop an- 
tenna functions [24] , where some attention has to be paid to the correct analytic continua- 
tion of the polylogarithmic functions. In the initial-final kinematics, the relevant final state 
phase space consists of three different Riemann sheets [53,54], which have to be patched 
together correctly. 

All quark- initiated initial-final antenna functions are listed in Table 1. We summarize 
their integrated forms in the following. 

5.1 Quark-quark antennae 

The next-to-leading order initial-final antenna functions were already integrated in [39] 
through to finite terms of order e^. Applying these antenna functions in an NNLO calcula- 
tion, one also needs their subleading terms through to e^. To this order, the quark- initiated 
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quark-antiquark antenna function is: 
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(5.1) 



At NNLO, the quark-initiated initial-final quark-antiquark double real radiation antenna 
are the crossings of the corresponding final-final quark-antiquark antenna functions. De- 
pending on the unresolved particles, they are of the A4-type (two unresolved gluons), B^- 
type (unresolved secondary quark pair of different flavour) and C4-type (pure interference 
contribution for secondary quark pair of identical flavour). 

The integrated form of the initial-final {q, ggq) antenna function reads at leading colour: 
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Crossing a primary quark (coupled to the external current) into the initial state, one 
obtains the {q, q'q'q) antenna function, whose integral yields 
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Correspondingly, crossing a secondary quark (not coupled to the external current) into 
the initial state, one obtains the {q' , qqq') antenna function, whose integral yields 
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+ 



r25 TT^ 



7 



20, 



Y - Y - 4C3 + y H(0; z) - -7r2H(0; z) + yH(l; z) - 8H(-1, 0; z) + H(0, 0; z) 

one 1 0Q (\A S. 

+ 2H(0, 1; z) + 4H(1, 0; z)+z-'{-^ + 2tt'' + -^H(0; z) + ^ H(l; z) - |h(-1, 0; z) 

- y H(0, 0; z) - y H(0, 1; z) - y H(l, 0; z) - y H(l, 1; z)) + 4H(1, 1; z) + 14H(0, 0, 0; z) 

+ 12H(0, 0, 1; z) + 8H(0, 1, 0; z) + 8H(0, 1, 1; z) + ^ (^ - ^ + ^H(l; z) 

8 , X 16 , s 16 , , 133 77r2 , ^ 

- -H(-l, 0; z) + y H(l, 0; z) + 1; z)) + z{-— + — - 4(3 - 4H(0; z) 

7 44 

- -7r2H(0; z) - — H(l; z) - 8H(-1, 0; z) - 11H(0, 0; z) - 14H(0, 1; z) - 4H(1, 0; z) 

- 4H(1, 1; z) + 14H(0, 0, 0; z) + 12H(0, 0, 1; z) + 8H(0, 1, 0; z) + 8H(0, 1, 1; z)) 
+ 0{e). 



(5.5) 



As explained in Section 2 above, there are three different crossings for the identical 
flavour only antenna functions C4, which are symmetrized over the antiquarks, but not 
over the quarks. 

Crossing either antiquark into the initial state yields the (g, qqq) antenna, whose inte- 
gral is: 



qmQ 



+- 



+ 



2 - ^ + ^H(0; .) - iH(0, 0; .) - 1h(1, 0; .) + ^ (^ 

+ H(1, 0;z))+z{~-^ + ^H(0; z) - iH(0, 0; z) - ^H(l, 0; z)) 
1 57r2 3C3 1_2 



+ -H(0;z) + H(0,0;2;) 



L4 12 



6 



7r^H(0; z) + 4H(1; z) + 2H(-1, 0; z) - 2H(0, 0; z) + H(0, 1; z) 



- 22 - 



- ^H(l, 0; z) - 2H(0, -1, 0; z) - ^H(0, 0, 0; z) - H(0, 0, 1; z) - H(0, 1, 0; z) 



7 



1 .n' 



- -H(l, 0, 0; z) - H(l, 0, 1; z) - H(l, 1, 0; z) + ^ (— + SCs + 3H(0; z) + -vr^HlO; z) 

+ ^H(0, 0; z) + ^H(0, 1; z) + 3H(1, 0; z) + 4H(0, -1, 0; z) + 3H(0, 0, 0; z) 
+ 2H(0, 0, 1; z) + 2H(0, 1, 0; z) + 7H(1, 0, 0; z) + 2H(1, 0, 1; z) + 2H(1, 1, 0; z)) 

+ + X^(0; ^) - ^^'H(0; z) - ^H(l; + 2H(-1, 0; z) + H(0, 0; z) 

+ H(0, 1; z) + ^H(l, 0; z) - 2H(0, -1, 0; z) - ^H(0, 0, 0; z) - H(0, 0, 1; z) - H(0, 1, 0; z) 



H(l, 0, 0; z) - H(l, 0, 1; z) - H(l, 1, 0; z) - ^) 



+ 0(e). 



(5.6) 



Crossing the quark which couples to the external current to the initial state yields the 
{q, qqq) antenna function, which integrates to: 



+- 



+ 



-5 + 2.= + Ml -)(-^ + ^ + ^ - ^) + f oW(^ - ^ + 2&) - 5H(0; .) 
- 6H(-1, 0; z) + 15H(0, 0; z) + 9H(1, 0; z) - ^H(-l, 0; -z) + ^ (-2 + 2H(0; z) 



1 



6H(-1, 0; z)) + z{27r^ - 2H(-1, 0; z) + 9H(0, 0; z) + 7H(1, 0; z)) + —— ( — 

6C3 + 10H(0; z) + ^7r2H(l; - 15H(0, 0; z) - 9H(1, 0; z) + 4H(0, -1, 0; z) 



+ 8H(l,0,0;z) +4H(l,l,0;z)) 
+ 0(6) . 



(5.7) 



Finally, crossing the quark not coupled to the external current (and thus participating 
in the coUinear splitting) , one obtains the (g, qqq) antenna function, integrating to: 



1 

+- 
e 



-2 - y - H(0; z) - 2H(-1, 0; z) + H(0, 0; z) + ^ (y + 4H(-1, 0; z) - 2H(0, 0; z)) 



+ z(2 + y - H(0; z) + 2H(-1, 0; z) - H(0, 0; z)) 

+ [ Y + y - 5C3 + 7r2H(-l; z) - ^H(0; z) - \t^^{^; z) - 4H(1; z) - 4H(-1, 0; z) 
- 2H(0, 0; z) - 2H(0, 1; z) + 4H(-1, -1, 0; z) - 8H(-1, 0, 0; z) - 4H(-1, 0, 1; z) 
-4H(0,-l,0;z) +3H(0,0,0;2) +2H(0,0,l;z) - ^H(-l,0;z) + ^(-2 + 2H(0;z) 



-23- 



+ 6H(-1, 0; z)) + ^ (+14C3 - ^7r2H(-l; z) - 4H(0; z) + ^7r^R{0; z) 
- 16H(-1, -1, 0; z) + 20H(-1, 0, 0; z) + 8H(-1, 0, 1; z) + 12H(0, -1, 0; z) 



11 TT^ 



19, 



- 6H(0, 0, 0; z) - 4H(0, 0, 1; z)) + z{-— + — + 5(3 - vr^Hl-l; z) + y H(0; z) 
+ j7r^H(0; z) + 4H(1; z) - 4H(-1, 0; z) - 2H(0, 0; z) - 2H(0, 1; z) - 4H(-1, -1, 0; z) 
+ 8H(-1, 0, 0; z) + 4H(-1, 0, 1; z) + 4H(0, -1, 0; z) - 3H(0, 0, 0; z) - 2H(0, 0, 1; z)) 
+ 0{e). (5.8) 



The integrated one-loop quark-antiquark antenna functions at leading and subleading 
colour read: 



1 



-^S{l-z) 



1 r 1 



53 



1 



+- 



L 6 



4-24^(^-^^-4 + 2^'^^^) 

13 w 43 57r\ 31 3 1 3H(0; ^) 

+ Sil~z){-- + —) + -Vo{z)-Viiz)--BiO;z)--ii{l;z)+ ^ > 



1 

+- 
e 



■z(^ + ^H(0;z) + ^H(l;z)) 
27 



16 24 
- 24 



2(1 -z) 



+ 



167 377r2 7C3x Stt^. , , 10_ , , _ , , 14^^, , 

-T + T + '^(^ - + ^ + 6 ) + - 12")^°^'^) - y ^^(^) + ^^("^^ - y ^(°' 

- ^H(l; z) - 2H(0, 0; 2f) - ^H(0, l-z)- ^H(l, 0; z) - H(l, 1; z) + ^ (^ 

49 5 25 TT^ 1 

+ -H(0; z) + 4H(0, 0; z) + -H(0, 1; z) + 3H(1, 0; z)) + z(- + — - -H(0; z) 

- ^H(l; z) - 2H(0, 0; z) - ^H(0, 1; z) - ^H(l, 0; z) - H(l, 1; z)) 
— r + 



67 677r2 29C3 w 773 237r2 eiTr^ 217(3 x / 



48 



48 



16 



1440 



36 



257 _ 137r2 



- ^)^o(.) + 
+ ^7r2H(l;z)- 



53 57r\ ^ ^ 29_ , , 2_ , , 117,,, , 
■y + — )^i(^) + -^^2{z) - -P3(^) - — H(0;z) 



12 

Vh(0;^) 



-H(l;z) 



119 , . 43 , 

— H(0,0;z)-yH(0,l;z) 



89 , . 71 , 
-H(l,0;z)--H(l,l;z) 



^H(0, 0, 0; z) - 3H(0, 0, 1; z) - 3H(0, 1, 0; z) - ^H(0, 1, 1; z) - 4H(1, 0, 0; z) 
^H(l, 0, 1; z) - 3H(1, 1, 0; z) + ^ (^ - ^ + |h(0; z) - ^K^^i^- z) 



+ i7r2H(l; 2) + ^H(0, 0; z) + ^H(0, 1; z) 

D 12 



19 



H(l,0;z) + 9H(0,0,0;z) 



+ 6H(0, 0, 1; z) + 6H(0, 1, 0; z) + -H(0, 1, 1; z) + 8H(1, 0, 0; z) + 5H(1, 0, 1; z) 



.1 57r2 29C3 49^ 



23, 



+ 6H(1, 1, 0; z)) + z(^ - ^ + ^ + gH(0; z) + -7r2H(0; z) + =^H(1; z) 



-24- 



+ ^7r2H(l; z) + ^H(0, 0; z) + ^H(0, 1; z) + 1h(1, 0; ^) + 1h(1, 1; z) 
- ^H(0, 0, 0; z) - 3H(0, 0, 1; z) - 3H(0, 1, 0; z) - ^H(0, 1, 1; z) - 4H(1, 0, 0; z) 
-^H(l,0,l;z) -3H(l,l,0;z) -2H(l,l,l;z)) -2H(1,1,1;^) 



(5.9) 



+- 
e 



1 + ^ + (5(1 - z) (-> ^) + - ^)Po(^) - 7H(0; z) - H(0, 0; 2:) - ^H(0, 1; z) 



5C; 



~8 6 



-H(0,0;z)-iH(0,l;z)) 



+ 



9C; 



61 257r2 437r'^ 



1440 + ^^^^) + (i - 5C3)^>o(^) + (-^ + \)V,{z) 



24 



35, 



16 96 

7, 



1 Q 7 1 

- — H(0; z) + — ^2jj(0; 2) - -H(l; z) + -7r2H(l; z) - 7H(0, 0; z) - -H(0, 1; z) 



24 

2H(1, 0; z) - |h(0, 0, 0; 2) - H(0, 0, 1; 2) - ^H(0, 1, 1; z) + ^H(l, 0, 0; z) 



+ iH(l,l,0;z) + ^| 



^8 2 4 ^'^ 12 



1 Q 1 

7r2H(0;z)--7r2H(l;z) 



+ 6H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z) + 5H(0, 0, 0; z) + 2H(0, 0, 1; 2) + ^H(0, 1, 1; z) 



-2H(l,0,0;z)-2H(l,l,0;z))+z(-J-^ + ^^ ' ^ 



13 

+ ^ + ^H(0;z) + — 7r2H(0;z) 



24 

+ ^H(l; z) + ^7r2H(l; z) - ^H(0, 0; z) + ^H(0, 1; z) - ^H(l, 0; z) - ^H(0, 0, 0; z) 
- H(0, 0, 1; z) - ^H(0, 1, 1; z) + ^H(l, 0, 0; z) + ^H(l, 1, 0; z)) 



+ 0{e) . 
The quark loop contribution is: 



(5.10) 



A, 



1 

<1,9Q 



+ ^[^^(1-^) 

+^i6 + 4^(l-)-3A)(.) + g 
+^ + - - ^) - ^^o(.) + ^P.(.) + ^H(0; .) + 1h(1; .) - ^ 

+ ^(-^ + ^H(0;z) + ^H(l;^)) 



H(0;2) 



-25- 



+ 



7_ 7r2 
6 16 



7C3^ 
9 ' 



+ (-i + ^)Vo{z) + \v,{z) ^ \v^{z) + 1h(0;z) 



12 12^ ' 4 ' ' 6 ' ' 2 
^H(l; z) + ^H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z) + ^H(l, 1; z) + 



6 



6 



6 



6 

1, 



1 



- ^H(0, 0; z) - ^H(0, 1; z) ^ ^H(l, 0; z)) + ^("^ " ^H(0; z) - ^H(l; z) + ^H(0, 0; z) 

+ iH(0, 1; z) + 0; + ^H(l, 1; z)) 

+ 0(6). (5.11) 
5.2 Quark-gluon antennae 

The integrated NLO quark-gluon initial-final antenna with an initial state quark reads 
through to for the gluon radiation: 

= 

1,99 



e 
1 

+ - 
e 



1 r 
1 



2(5(1 - z) 



11 



+ 



+e 



1 + —5{1-z) + z-2Vq{z) 

1 + 5(1 - z) (^ - y ) - ^ - ^Po(^) + 2Pi {z) + H(0; + H(l; z) - y4^H(0; 

+ z(-l + H(0;2) +H(l;z)) 
^ vr^ ^ ,202 Il7r2 MCsn . 67 

+ H(l; z) + H(0, 0; z) + H(0, 1; z) + H(l, 0; ^) + H(l, 1; ^) - 1 (^ + iH(0; z) + z)) 

+ (-^H(0; z) - 2H(0, 0; z) - 2H(0, 1; - 2H(1, 0; z)) + ^-^ - H(0; 
1 — z 6 ' ^ 4 

- H(l; + H(0, 0; z) + H(0, 1; z) + H(l, 0; z) + H(l, 1; z)) 

7C3^.,^ wl214 677r2 77^,^ 

T-^+'^^'-")(^-^-48 

67 7r\„ . . 11„ _ 1_ . . . 1 



202 llvr^ 14C: 



18^+^°(^)(-^+ 24 



+ 



4 



+ ( 18 ~ V)^i(^) - T^^2(^) + iD^{z) + H(0; z) - -7r2H(0; z) + H(l; z) - -^Hil; z) 



12 



+ H(0, 0; z) + H(0, 1; z) + H(l, 0; z) + H(l, 1; z) + H(0, 0, 0; z) + H(0, 0, 1; z) 

+ H(0, 1, 0; z) + H(0, 1, 1; z) + H(l, 0, 0; z) + H(l, 0, 1; z) + H(l, 1, 0; z) + H(l, 1, 1; z) 



1 



13 TT^ 



1. 



1. 



1, 



; 27 + 12 " 9^^°' " 9^^^' " 3^^°' °' " 3^^°' ^' " 3^^^' °' 

- iH(l, 1; z)) + ^ (-^H(0; z) + i7r2H(0; z) - Hh(0, 0; z) - Hh(0, 1; 

- — H(l, 0; z) - 2H(0, 0, 0; z) - 2H(0, 0, 1; z) - 2H(0, 1, 0; z) - 2H(0, 1, 1; z) 
6 

- 2H(1, 0, 0; z) - 2H(1, 0, 1; z) - 2H(1, 1, 0; z)) + z(^ - ^ - ^7r2H(0; z) 

- ^vr^H(l; z) - H(0, 0; z) - H(0, 1; z) - H(l, 0; z) - H(l, 1; z) + H(0, 0, 0; z) 



-26- 



+ H(0, 0, 1; z) + H(0, 1, 0; z) + H(0, 1, 1; z) + H(l, 0, 0; z) + H(l, 0, 1; z) + H(l, 1, 0; z) 
+ H(l,l,l;z)) 

+ 0{e^), (5.12) 

while we have to distinguish two cases for the secondary quark radiation. Crossing the 
primary quark into the initial state, one obtains 



1 

e 



-5{l-z) 



+e 



164 Stt^ TCsx /28 tt' 



H(0;z) 



V[+^(l - .)(-- + ^ + f ) + - T^)^^o(.) - -^V,iz) + + ^(gHlO;.) 

1 1 1 1 26 vr^ 4 4 

+ -H(0, 0; + -H(0, 1; z) + -H(l, 0; z)) +-(-- + - - -H(0; z) - -H(l; z) 



- ^H(0, 0; z) - ^H(0, 1; ^) - ^H(l, 0; z) - ^H(l, 1; z)) 
while crossing the secondary quark results in 



(5.13) 



cO _ 

^q',q'q " 

Ir ^ 1 Z 
+ - -1 + - + - 

e L z 1 

+ [-\- H(0; z) - H(l; z) + -^{2 + H(0; z) + H(l; z)) + z(^H(0; z) + ^H(l; z)) 

+e [-3 + ^ - ^H(0; z) - ^H(l; z) - H(0, 0; z) - H(0, 1; z) - H(l, 0; z) - H(l, 1; z) 
1 7r2 

+ - (4 - — + 2H(0; z) + 2H(1; 2) + H(0, 0; z) + H(0, 1; 2) + H(l, 0; z) + H(l, 1; 2)) 
+ ^(-y + ^H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z) + ^H(l, 1; z)) 
+^ [-6 + ^ + ^ - 3H(0; ^) + ^7r2H(0; ^) - 3H(1; z) + ^7r2H(l; ^) - ^H(0, 0; z) 

- ^H(0, 1; z) - ^H(l, 0; z) - ^H(l, 1; z) - H(0, 0, 0; z) - H(0, 0, 1; z) - H(0, 1, 0; z) 

1 TT^ 

- H(0, 1, 1; z) - H(l, 0, 0; z) - H(l, 0, 1; z) - H(l, 1, 0; z) - H(l, 1, 1; z) + - (8 - y 

_ + 4H(0; z) - ^7r2H(0; z) + 4H(1; z) - ^vr2H(l; z) + 2H(0, 0; z) + 2H(0, 1; ^) 
+ 2H(1, 0; z) + 2H(1, 1; z) + H(0, 0, 0; z) + H(0, 0, 1; z) + H(0, 1, 0; z) + H(0, 1, 1; z) 



-27- 



+ H(l, 0, 0; z) + H(l, 0, 1; z) + H(l, 1, 0; z) + H(l, 1, 1; z)) + z{--^ - j7r2H(0; z) 



- ^7r2H(l; z) + iH(0, 0, 0; z) + hl{Q, 0, 1; z) + ^H(0, 1, 0; z) + 1, 1; z) 

+ ^H(l, 0, 0; z) + ^H(l, 0, 1; ^) + ^H(l, 1, 0; z) + ^H(l, 1, 1; z)) 

+ O {e^) . (5.14) 

The NNLO double real radiation quark-gluon antenna functions are of the type 
(quark-gluon-gluon-gluon) and of £?4-type (quark-quark-antiquark-gluon). While the for- 
mer has only one possible quark-initiated crossing, we have to distinguish several different 

cases for the latter. 

The integrated form of the initial-final (g, ggg) antenna function reads: 

= 

1,999 



1 r 
1 

+ -7 



+ - 



+- 



— z 



4 

z 



10(5(1 - z) 

10 + 22(5(1 -z)- 2QVo{z) + lOz 

17 + (5(1 -z){^- Stt^) - UVq{z) + 40X»i(z) + 22H(0; z) + 20H(1; - 

9 1 

+ z(5 + 22H(0; z) + 20H(1; z)) 

107 257r2 ,,4868 1547r2 272C3x , 1180 ^ , , , , / x 

^ ^ + (5(1 - z) {— f-) + ( — + 167r2)Po(^) + 88Pi(z) - 40P2(^) 

113 

H(0; z) + 34H(1; z) + 50H(0, 0; z) + 44H(0, 1; z) + 42H(1, 0; z) + 40H(1, 1; z) 

3 

+ ^ (-¥ - ^H(0; - 88H(0, 0; z) - 80H(0, 1; z) - 84H(1, 0; z)) + -(-^ 



1-z 



8H(0;z) -8H(l;z)) + z{22 



257r2 41 



+ — H(0; + 10H(1; z) + 50H(0, 0; z) 



+ 44H(0, 1; z) + 42H(1, 0; z) + 40H(1, 1; z)) 



+ 



565 1457r2 3O8C3 



L 9 
+ 



18 
.2 



3 ^ ^^54 
,2360 



25811 14087r2 31971^ 934C; 



27 



+ 



180 
80, 

y 



-)+2?o(^)(- 



9736 



27 



616, 



50 2 



7r^H(l; z) - 8H(-1, 0; z) + 121H(0, 0; z) 



3 ^ ^ 9 

-|vr2H(0;.) + ^H(l;. 
226 

— H(0, 1; z) + 96H(1, 0; .z) + 68H(1, 1; z) - 16H(-1, -1, 0; z) + 8H(-1, 0, 0; z) 
+ 4H(0, -1, 0; z) + 106H(0, 0, 0; z) + 100H(0, 0, 1; z) + 84H(0, 1, 0; z) + 88H(0, 1, 1; 2) 

+ 82H(l,0,0;z) +84H(l,0,l;z) +80H(l,l,0;z) +80H(l,l,l;z) + ^(-^ + ^ 

2 224 224 2 

- 4C3 - 37r2H(-l; 2) - — H(0; z) - — H(l; z) + -7r2H(l; - 4H(-1, 0; z) 

44 

- 16H(0, 0; z) - 16H(0, 1; z) - — H(l, 0; z) - 16H(1, 1; z) - 8H(-1, -1, 0; z) 



-28- 



+ 4H(-1, 0, 0; z) + 8H(0, -1, 0; z) + 8H(1, 0, 0; z) + 4H(1, 1, 0; z)) + - 
2368„,„ , 116 2TT/n N 4 



llTT^ 



32C3 - -^H(0; z) + ^7r^H(0; z) + |7r^H(l; z) - 243H(0, 0; z) - ^H(0, 1; z) 



z 

554, 
3" 



620, 



- ^H(l, 0; z) + 8H(0, -1, 0; z) - 184H(0, 0, 0; z) - 176H(0, 0, 1; z) - 152H(0, 1, 0; z) 

443 

- 160H(0, 1, 1; z) - 164H(1, 0, 0; z) - 168H(1, 0, 1; z) - 160H(1, 1, 0; z)) + z{— 



257r^ 

18 

50 



320C3 



391 



65 



-7r^H(-l; z) + — H(0; z) - — 7r^H(0; z) + 44H(1; z) 



82, 



— 7r^H(l; z) - 4H(-1, 0; z) + 47H(0, 0; z) + — H(0, 1; z) + 34H(1, 0; z) 
3 3 

+ 20H(1, 1; z) - 8H(-1, -1, 0; z) + 4H(-1, 0, 0; z) + 106H(0, 0, 0; z) + 100H(0, 0, 1; z) 
+ 84H(0, 1, 0; z) + 88H(0, 1, 1; + 82H(1, 0, 0; z) + 84H(1, 0, 1; z) + 80H(1, 1, 0; z) 

+ 80H(l,l,l;z)) 

+ 0(e). (5.15) 



For the E^-type antenna functions, we can either cross the primary quark or the 
secondary quark into the initial state. In both cases, we have leading and subleading 
colour contributions. 

Crossing the primary quark, one obtains the {q, q'q'g) antenna functions, whose integral 
reads at leading colour 



q,q'q'9 



1 



47 10 2 ■ 

-1 + 



If 4 X/ 55 7 237r2 94_ , , 20_ , , 7,,, , , 22H(0;z) 

1 26 4 4 3 7 

+ -(--- -H(0; z) - -H(l; z))+z{--- -H(0; z) - 2H(1; z)) 

- + -+5(l-.)(-^ + ^+ 9 - M8 9 



+ 



28613 , 257r2 , 118C3^ ^ .557 237r^. 188 . . ^ 20 

) + {-^ - -^Po{z) - —Vi{z) + —V2{z) 



38 8 19 14 16 

- -H(0; z) - -H(l; z) - -H(0, 0; z) - -H(0, 1; z) - -H(l, 0; z) - 4H(1, 1; z) 

1 , , , 50 , . 44 , . 44 , 1 , 103 57r2 

+ ^ (22H(0; z) + y H(0, 0; z) + y H(0, 1; z) + y H(l, 0; z)) +-(-— + — 

52 52 8 8 8 8 

- -H(0; z) - -H(l; z) - -H(0, 0; z) - -H(0, 1; z) - -H(l, 0; z) - -H(l, 1; z)) 

+ ^(-T^ + ^ - f H(0; z) - 3H(1; z) - ^H(0, 0; z) - ^H(0, l;z)- ^H(l, 0; z) 



12 

-4H(l,l;z)) 



(5.16) 



-29- 



and at subleading colour 



+- 



+- 



+- 



+ 



cO 

^q,q'q'g 

1 



-3*(1-.) 

-f*(l-.) + iPo(.)-£^ 

38 



373 Svr^ 



s , , 8^ , , 8H(0;z) 



1 / 25 4 , , 4 , 



6973 957r2 64C3x /373 IOtt^. , , 76_ , , 8_ , , 1 /76,,,^ , 
+ f H(0, 0; .) + f H(0, 1; .) + f H(l, 0; .)) + ^ (-^ + ^ 

en en 8 8 8 8 

- -H(0; z) - -H(l; z) - -H(0, 0; z) - -H(0, 1; z) - -H(l, 0; z) - -H(l, 1; z)) 



+ 0(6). 



(5.17) 



Crossing the secondary quark, one obtains the (q', q'qg) antenna functions, whose in- 
tegral reads at leading colour 



cO 

q',q'Q9 



+ 



+ 



-4+- + 2Z 
z J 

42;2 1 

15 - — - 4H(0; z) - 8H(1; z) + -{— + 8H(0; z) + 8H(1; z)) 

Z o 



3 

+ ^(^ + fH(0;z)+4H(l;^)) 
1 r 23 5 77r^ 

+ — + ^ - 39H(0; z) - 30H(1; z) + 4H(-1, 0; z) - 8H(0, 0; z) - 8H(0, 1; z) 

- 16H(1, 0; z) - 16H(1, 1; z) + - - Stt^ + 28H(0; z) + ^H(l; z) + 4H(-1, 0; z) 

47 77r^ 5 

+ 16H(0, 0; z) + 16H(0, 1; z) + 16H(1, 0; z) + 16H(1, 1; z)) + z{— - — - -H(0; z) 

X ^ o ^ 

31 

+ 3H(1; z) + 2H(-1, 0; z) + yH(0, 0; z) + 13H(0, 1; z) + 8H(1, 0; z) + 8H(1, 1; z)) 

+ -yH(0;.)-^H(l;z)) 

r 2695 1157r2 17OC3 o / x 337 , , 14 2 . ^ 235 , , 
+ [-^ + + - 27r2H(-l; z) - — H(0; z) + y tt^HCO; z) - — H(l; z) 

20 

+ — 7r2H(l; z) - 89H(0, 0; z) - 78H(0, 1; z) - 52H(1, 0; z) - 60H(1, 1; z) 

- 8H(-1, -1, 0; z) + 16H(-1, 0, 0; z) + 8H(-1, 0, 1; z) + 8H(0, -1, 0; z) - 16H(0, 0, 0; z) 

- 16H(0, 0, 1; z) - 12H(0, 1, 0; z) - 16H(0, 1, 1; z) - 32H(1, 0, 0; z) - 32H(1, 0, 1; z) 



- 32H(1, 1, 0; z) - 32H(1, l,l;z) + -{ 
+ 1^H(0; z) - 8n'm z) + ^H(l; z) - 



1 .37925 9l7r2 134(3 



6 



27r^H(-l;z) 



y7r2H(l; z) + ^H(-l, 0; + 56H(0, 0; z) 



-30- 



205 220 

+ 56H(0, 1; z) + — H(l, 0; z) + — H(l, 1; z) - 8H(-1, -1, 0; z) + 16H(-1, 0, 0; z) 

+ 8H(-1, 0, 1; z) + 8H(0, -1, 0; z) + 32H(0, 0, 0; z) + 32H(0, 0, 1; z) + 28H(0, 1, 0; z) 
+ 32H(0, 1, 1; z) + 32H(1, 0, 0; z) + 32H(1, 0, 1; z) + 32H(1, 1, 0; z) + 32H(1, 1, 1; z)) 



1069 37r2 82(3 



83 



47, 



^72 4 ~ ^^(-^-^ ^) + 22H(0; z) - ^7r^H(0; z) + -H(l; z) 



10 



23 



— 7r^H(l; z) - 4H(-1, 0; z) - — H(0, 0; z) - 5H(0, 1; z) + 3H(1, 0; z) + 6H(1, 1; z) 

67 

- 4H(-1, -1, 0; z) + 8H(-1, 0, 0; z) + 4H(-1, 0, 1; z) + 4H(0, -1, 0; z) + yH(0, 0, 0; z) 
+ 31H(0, 0, 1; z) + 24H(0, 1, 0; z) + 26H(0, 1, 1; z) + 16H(1, 0, 0; z) + 16H(1, 0, 1; z) 
+ 16H(1, 1, 0; z) + 16H(1, 1, 1; z)) +z''{-— + 27r^ + — H(0; z) + -H(l; ^) 



2H(l;z)) 



- ^H(-l, 0; z) - y H(0, 0; z) - ^ H(0, 1; z) - ^H(l, 0; z) - y H(l, 1; z)) 
+ 0{e) 

and at subleading colour 

1 r 2 

+^ -2 + - + Z 

+4 f-8 - 5H(0; z) - 4H(1; z) + - (9 + 4H(0; + 4H(1; z)) + 2(7 + ^H(0; , 
1 r 127 

+- [ — ^ + 27r2 - 16H(0; z) - 16H(1; z) - 11H(0, 0; z) - 10H(0, 1; z) - 8H(1, 0; z) 
1 75 Stt^ 

- 8H(1, 1; -z) + - (17 - ^ + 18H(0; z) + 18H(1; + 8H(0, 0; z) + 8H(0, 1; z) 

7 19 7 11 

+ 8H(1, 0; z) + 8H(1, 1; z)) + z(- - tt^ + -H(0; z) + -H(l; z) + y H(0, 0; z) 

+ 5H(0, 1; z) + 4H(1, 0; z) + 4H(1, 1; z)) 

r 499 207r2 7OC3 259 , , 9 , , , 127 , , 10 o , , 
+ [-— + ^ + ^- — H(0; z) + -TT^HIO; z) - — H(l; z) + yTr^HCl; z) 

- 32H(0, 0; z) - 32H(0, 1; z) - 32H(1, 0; z) - 32H(1, 1; z) - 23H(0, 0, 0; z) 

- 22H(0, 0, 1; z) - 20H(0, 1, 0; z) - 20H(0, 1, 1; z) - 16H(1, 0, 0; z) - 16H(1, 0, 1; z) 



(5.18) 



- 16H(1, 1, 0; z) - 16H(1, 1, 1; z) + 
10 



1 ,565 157r2 „ ^ 10 o.w x 

-(^ ^ + 75H(0;z)-y7r2H(0;z) 



+ 75H(1; 2) - — 7r^H(l; 2) + 36H(0, 0; z) + 36H(0, 1; 2) + 36H(1, 0; z) + 36H(1, 1; z) 
+ 16H(0, 0, 0; z) + 16H(0, 0, 1; z) + 16H(0, 1, 0; z) + 16H(0, 1, 1; z) + 16H(1, 0, 0; z) 



+ 16H(1, 0, 1; z) + 16H(1, 1, 0; z) + 16H(1, 1, 1; z) 



+ 7H(0; z) 



- ^7r2H(0; 2) + 7H(1; z) - ^tt^HII; z) + ^H(0, 0; z) + y H(0, 1; z) + 7H(1, 0; z) 
23 

+ 7H(1, 1; z) + yH(0, 0, 0; z) + 11H(0, 0, 1; z) + 10H(0, 1, 0; z) + 10H(0, 1, 1; z) 



-31 - 



+ 8H(1, 0, 0; z) + 8H(1, 0, 1; z) + 8H(1, 1, 0; z) + 8H(1, 1, 1; z) 
+ 0{e). 



35C3 



(5.19) 



The one-loop virtual antenna functions at NNLO are obtained by crossing the results 
from [24]. For the {q, gg)-ca.se, we have the leading colour and the quark loop correction, 
which integrate to 



q,99 

1 



1 



-2^(1-) 



' ^ ^ + ^(1 (-1 + ^) + y ^o(.) - 2V, iz) - 2H(0; .) - H(l; .) + + j- 



+ z(---2H(0;z)-H(l;z)) 



1 

+- 
e 



17 77r2 ^, 



1967 36l7r2 



108 



+ 



72 



22, 

y 



23 17 

- yH(0; z) - y H(l; z) - 6H(0, 0; z) - 3H(0, 1; z) - 3H(1, 0; z) - 2H(1, 1; z) 

1 1 14 4 

+ (11H(0; z) + 12H(0, 0; z) + 6H(0, 1; 2) + 6H(1, 0; z)) + - (-5- + oH(0; z) 

+ ^H(l; z))+z{^+^-^ + iH(0; ^) - ^H(l; z) - 6H(0, 0; z) - 3H(0, 1; z) 



+ 



-3H(1,0;^) -2H(l,l;z)) 
17 677r2 25C3 w 



14453 30237r2 
+ 



324 



432 



llvr-^ 269C3X /3197 527r2 
+ — ^) + (- 



72 



6 



108 



9 



^)V,{z) + (-^ + ^)Pi(z) + yP2(^) - \v,{z) - ^H(0;z) + ^7r2H(0;z) - ^H(l;z) 



4 „ 79 29 

+ -7r2H(l; z) - ^H(0, 0; z) - -|h(0, 1; z) 



^H(l,0;^)-fH(l,l;2) 



- 14H(0, 0, 0; z) - 8H(0, 0, 1; z) - 6H(0, 1, 0; z) - 5H(0, 1, 1; z) - 6H(1, 0, 0; z) 

- 5H(1, 0, 1; z) - 4H(1, 1, 0; z) - 4H(1, 1, 1; z) + ^ _ 2C3 + ^H(0; z) 

- y TT^HCO; z) - ^TT^HCl; z) + y H(0, 0; z) + y H(0, 1; z) + y H(l, 0; 2;) 

+ 28H(0, 0, 0; z) + 16H(0, 0, 1; z) + 12H(0, 1, 0; z) + 10H(0, 1, 1; z) + 12H(1, 0, 0; z) 

1 145 llvr^ 53 5 
+ 10H(1, 0, 1; z) + 8H(1, 1, 0; z)) + - (— - ^ + -H(0; z) + -H(l; z) 

+ y H(0, 0; + 2H(0, 1; z) + -H(l, 0; z) + -H(l, 1; z)) + z(y^ + ^ + y ^(0; z) 
+ ^7r2H(0; z) + ^ H(l; z) + ^7r2H(l; z) + ^ H(0, 0; z) + ^H(0, 1; z) - ^H(l, 0; z) 

+ ^H(l, 1; z) - 14H(0, 0, 0; z) - 8H(0, 0, 1; z) - 6H(0, 1, 0; z) - 5H(0, 1, 1; z) 



-32- 



- 6H(1, 0, 0; z) - 5H(1, 0, 1; z) - 4H(1, 1, 0; z) - 4H(1, 1, 1; z)) 
+ 0{e) 

and 



(5.20) 



= 

<1,99 

+ z(-^ + ^H(0;z) + ^H(l;z)) 



+ - 



1 



z' 
+ 3 



1 

9z 



+ 



vr 



. ,,979 llvr^ MCsx , 38 7r\ , x -../^ x 

2 - T + (3M - ^ - ^) + |- 27 + t)^;(^) + I8^^(^) - 3^^^J^) + 3«(«' 

+ ^H(l; z) + ^H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z) + ^H(l, 1; z) + ^ (-^H(0; z] 

- ^H(0, 0; z) - ^H(0, 1; z) - ^H(l, 0; z)) + 1 (- A - 1h(0; - 1h(1; z)) + 

- iH(0; z) - ^H(l; z) + ^H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z) + ^H(l, 1; z)) 



(5.21) 



For the {q,q'q')-case (primary quark in the initial state), there are leading and sub- 
leading colour as well as quark loop contributions. Their integrals are: 



1 



+ - 



+ 



1 rii9 

108 

1 



5(1 -z) 



11 , , 2H(0;z) 1 ,2 1 
18^°(^)-3(rrin(9 + 3^(°=^)) 



6 + - ^Hil - ^) + ^ - i^o(.) + ^P.(.) - ^H(0; z) + ^H(0, 0; z) 



+ -H(l,0;z) + - 



.(_^_-H(0;, 
^9 18 ^ ' 



^H(0,0;z)-^H(0,l;z)-^H(l,0;z)) 



+ ^ - ^ + ^H(0; z) - 1h(1; z) + H(0, 0; z) + 1; ^) + ^H(l, 0; z)) 



+ -2(^ + ^H(0,0;z) + iH(l,0;z)) 
+ 0(e), 



(5.22) 



+ 



1 rl 
?l-6 



5(1 -z) 



-33- 



1 

1 

+- 
e 



|*(i-.)-ii,„(.)+i; 



173 57r^ 



19, 



+ 



- ^H(0, 0; ^) - ^H(0, 1; z) - ^H(l, 0; ^)) + i - + ^R{0; z) + |h(1; ^) 
+ ^H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z) + ^H(l, 1; z)) 



173 57r2, 



19, 



19, 



+ 0{e) 
and 



(5.23) 



+ - 



1 r 5 



+ 



,,Sil-z)--Voiz) + - 
,53 



1 



H(0;2) 



1/1 1, 



1, 



+ 0(6). 



(5.24) 



Likewise, the {q' , q' q)-case (secondary quark in the initial state), there are leading and 
subleading colour as well as quark loop contributions, which integrate to: 

- 



1 

+ - 
e 



- + H(0; z) + -H(l; z) + -{--- H(0; z) - -H(l; z)) - z{- + -H(0; z) + -H(l; z)) 

41 TT^ 5 7 3 

- y + 3H(0; z) + -H(l; z) + 4H(0, 0; z) + 2H(0, 1; z) + 2H(1, 0; z) + -H(l, 1; z) 

1 29 TT^ 8 9 
+ - (gg + y - 3H(0; ^) - -H(l; z) - 4H(0, 0; z) - 2H(0, 1; z) - 2H(1, 0; ^) 

- ^H(l, l;z))+z{'^ + ^ + ^H(0; z) - ^H(l; z) - 2H(0, 0; z) - H(0, 1; z) 

-H(l,0;z)-^H(l,l;z)) 



+ 



347 437r2 



5C3 - :^H(0; z) - ^vr2H(0; z) + ^H(l; z) - ^7r2H(l; z) + ^H(0, 0; z) 
6b b 3d 4 z 



54 12 

+ ^H(0, 1; z) - ^H(l, 0; z) + ^H(l, 1; z) + 10H(0, 0, 0; z) + 6H(0, 0, 1; z) 

O O t: 

+ 4H(0, 1, 0; z) + 4H(0, 1, 1; z) + 5H(1, 0, 0; z) + 4H(1, 0, 1; z) + 4H(1, 1, 0; z) 

3 9 . . 15 , X 17 , . 19 , X 21 , 
+ ^7r2H(l; z) - y H(0, 0; z) - y H(0, 1; ^) - y H(l, 0; z) - y H(l, 1; z) 

- 10H(0, 0, 0; z) - 6H(0, 0, 1; z) - 4H(0, 1, 0; z) - 4H(0, 1, 1; z) - 5H(1, 0, 0; z) 

- 4H(1, 0, 1; z) - 4H(1, 1, 0; z) - ^H(l, 1, 1; ^)) + ^(^ + ^ + ^ + ^H(0; z) 



-34- 



+ ^7r2H(0; z) + ^H(l; z) + ^7r2H(l; z) + ^H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z) 

- |h(1, 1; z) - 5H(0, 0, 0; z) - 3H(0, 0, 1; z) - 2H(0, 1, 0; z) - 2H(0, 1, 1; z) 
8 

-^H(l,0,0;z) -2H(l,0,l;z) -2H(l,l,0;z) -^H(l,l,l;^)) 



+ 0(e), 



(5.25) 



1 _ 1 _ z 

2 ~ 2i ~ 4 
7 



1 

1 



^ + H(0; + 1h(1; ^) " ^ + H(0; z) + z)) - z{^ + z) + ^H(l; z)) 



+ - 
e 



1 rl3 TT^ 



7, 



Y - Y + 2H(0; 2;) + -H(l; z) + 2H(0, 0; z) + H(0, 1; z) + H(l, 0; z) + -H(l, 1; z) 



1, 35 TT^ 9^,, 

+ ^(-T + T-2«^o^' 



-H(l; - 2H(0, 0; z) - H(0, 1; z) - H(l, 0; 



-iH(l,l;z))+<-l + Y-^H(0;z) 



^H(l;2;)-H(0,0;z)-^H(0,1;2;) 



-^H(l,0;z)-^H(l,l;z)) 



+ 



^H(l,0;z)-i 

43 Tvr^ 13(3 
T 

7„.„ . . 7 



^ + 13H(0; z) - ^TT^HCO; z) + y H(l; z) - ^7r2H(l; z) + 7H(0, 0; 

7, 



+ -H(0, 1; z) + -H(l, 0; z) + -H(l, 1; z) + 4H(0, 0, 0; z) + 2H(0, 0, 1; z) + 2H(0, 1, 0; z) 
+ H(0, 1, 1; z) + 2H(1, 0, 0; z) + H(l, 0, 1; + H(l, 1, 0; z) + ^H(l, 1, 1; z) + ^ (-^ 



+ ?|!^I|3_|H(0;.) + i.^H(0;.) 



35 1 

— H(l;z) + -7r2H(l;z)-9H(0,0;z) 



- ^H(0, 1; z) - ^H(l, 0; z) - ^H(l, 1; z) - 4H(0, 0, 0; z) - 2H(0, 0, 1; z) - 2H(0, 1, 0; z) 

- H(0, 1, 1; z) - 2H(1, 0, 0; z) - H(l, 0, 1; z) - H(l, 1, 0; z) - ^H(l, 1, 1; z)) + z{-2 
+ ^ + ^ - 2H(0; z) + ^TT^BiO; z) - H(l; z) + ^tt^HII; z) - ^H(0, 0; z) 

- ^H(0, 1; z) - ^H(l, 0; z) - ^H(l, 1; z) - 2H(0, 0, 0; z) - H(0, 0, 1; z) - H(0, 1, 0; z) 
-^H(0,l,l;z) -H(1,0,0;^) -^H(1,0,1;^) -^H(l,l,0;z) -^H(l,l,l;^)) 



+ 0{e) 
and 



(5.26) 



+- 



+ 



-35- 



^ - ^H(0; z) - ^H(l; z) - H(0, 0; z) - hl{Q, 1; z) - 1h(1, 0; z)) 



14 TT^ 



1, 



1. 



1, 



+ 0(6). 



(5.27) 



5.3 Gluon-gluon antennae 

The gluon-gluon antenna functions can contribute to quark-initiated processes, since they 
also contain configurations with the splitting of a gluon into a quark-antiquark pair. 
At NLO, one has only the (g, qg) gluon-gluon antenna function, whose integral is 



1 
e 

+ 
+e 



-■ 1 -2 
■' + i + 2 

-1 - H(0; z) - H(l;2) + i( J + H(0; z) + H(l; 2)) + z(\-a(0; z) + ^HCl; z)) 



vr 



-2 + — - H(0; z) - H(l; z) - H(0, 0; z) - H(0, 1; z) - H(l, 0; z) - H(l, 1; z) 
+ ; - ^ + Jh(0; z) + ^H(l; z) + H(0, 0; z) + H(0, 1; z) + H(l, 0; z) + H(l, 1; z)) 
+ ^(-y + ^H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z) + ^H(l, 1; 2)) 

-4 + ^ + ^ - 2H(0; z) + ^7r2H(0; z) - 2H(1; ^) + ^vr^HCl; z) - H(0, 0; z) 

- H(0, 1; 2) - H(l, 0; z) - H(l, 1; - H(0, 0, 0; z) - H(0, 0, 1; z) - H(0, 1, 0; z) 

- H(0, 1, 1; z) - H(l, 0, 0; z) - H(l, 0, 1; z) - H(l, 1, 0; z) - H(l, 1, 1; 2;) + i (7 - ^ 

- ^ + ^H(0; - ^7r2H(0; z) + ^H(l; z) - ^tt^HCI; + ^H(0, 0; z) + ^H(0, 1; z) 
+ ^H(l, 0; z) + ^H(l, 1; z) + H(0, 0, 0; z) + H(0, 0, 1; z) + H(0, 1, 0; z) + H(0, 1, 1; z) 
+ H(l, 0, 0; z) + H(l, 0, 1; + H(l, 1, 0; z) + H(l, 1, 1; z)) + - ^7r2H(0; z) 

DO 



1 



1 



1, 



1, 



- g7r^H(l; z) + -H(0, 0, 0; z) + -H(0, 0, 1; z) + -H(0, 1, 0; z) + -H(0, 1, 1; z 

+ ^H(l, 0, 0; z) + ^H(l, 0, 1; z) + ^H(l, 1, 0; z) + ^H(l, 1, 1; z)) 
+ 0(e3) . 



(5.28) 



At NNLO, the double real radiation G^-type and i?4=type antenna functions allow 
for crossings with initial state quarks. 

The (g, qgg) at leading colour integrates to 



-.0 

^1,199 



+ 



-8+- + 4z 
z 



-36- 



1 r 82 



1 ,104 



,53 21 



_ _ 9H(0;z) - 16H(l;z) + -(^ + 16H(0;z) + 16H(1; z)) + z(- + -H(0;2:) 



+ 8H(l;z)) 



1 

+- 
e 



+ 



+ ^ - ^H(0; z) - ^H(l; z) + 4H(-1, 0; z) - 15H(0, 0; z) 

- 18H(0, 1; z) - 32H(1, 0; z) - 32H(1, 1; z) + -(^ - Gtt^ + ^H(0; z) + ^H(l; z) 

z b 6 6 

301 

+ 4H(-1, 0; z) + 32H(0, 0; z) + 32H(0, 1; z) + 32H(1, 0; 2) + 32H(1, 1; z)) + z{— 

2377-2 c-o 

- + Y^H(0; z) + -H(l; z) + 2H(-1, 0; z) + y H(0, 0; z) + 21H(0, 1; z) 

00 IK 8 

+ 16H(1, 0; z) + 16H(1, 1; z)) y H(0; z) - -H(l; z)) 

6625 2567r2 286(3 „ 2... x 3181 , 55 n,,, , 
- - 27r^H(-l; - — — H(0; + — 7r2H(0; z) 



27 



+ 



9 



+ 



36 



6 



2291 

18 



H(l;z) 



oQ S73 
+ yTT^Hll; z) - 8H(-1, 0; z) - — H(0, 0; ^) - :^H(0, 1; ^) - ^H(l, 0; z) 

328 

- — H(l, 1; z) - 8H(-1, -1, 0; z) + 16H(-1, 0, 0; z) + 8H(-1, 0, 1; z) + 16H(0, -1, 0; z) 

- 27H(0, 0, 0; z) - 30H(0, 0, 1; z) - 28H(0, 1, 0; z) - 36H(0, 1, 1; z) - 60H(1, 0, 0; z) 

- 64H(1, 0, 1; .) - 64H(1, 1, 0; .) - 64H(1, 1, 1; .) + i (-^^^^^ 



370, 

T 



298, 

T 



z^ 108 



38 



- 27r2H(-l; z) + ^H(0; 2) - leTr^HCO; z) + ^H(l; - y^r^HCl; z) 

- |h(-1, 0; .) + ^H(0, 0; .) + ^H(0, 1; .) + ^H(l, 0; z) + ^H(l, 1; z) 

- 8H(-1, -1, 0; z) + 16H(-1, 0, 0; z) + 8H(-1, 0, 1; z) + 64H(0, 0, 0; z) + 64H(0, 0, 1; z) 
+ 56H(0, 1, 0; z) + 64H(0, 1, 1; z) + 60H(1, 0, 0; z) + 64H(1, 0, 1; z) + 64H(1, 1, 0; z) 



+ 64H(1,1,1;2)) +z 
131 . 301 

Is 



7411 1997r2 149C3 Ott/ . x 511 , 
^-vr'H(-l;2;) + — H(0;2;) 



216 

19 65 
7r2H(0; + — H(l; z) - —^^^{1; z) + — H(0, 0; z) + ^H(0, 1; z) 



72 
19 

y 



61. 
"6" 



12 

50 53 

+ y H(l, 0; z) + y H(l, 1; 2) - 4H(-1, -1, 0; z) + 8H(-1, 0, 0; z) + 4H(-1, 0, 1; z) 
99 

+ y H(0, 0, 0; z) + 47H(0, 0, 1; z) + 38H(0, 1, 0; z) + 42H(0, 1, 1; z) + 30H(1, 0, 0; z) 

208 1 28 

+ 32H(l,0,l;z) +32H(l,l,0;z) +32H(l,l,l;z)) + 27r2 + — H(0;z) 

+ ^H(l; - ^H(-l, 0; z) - ^H(0, 0; z) - f H(0, 1; z) - ^H(l, 0; z) - ^H(l, 1; z)) 



9 

+ 0{e) , 
while the subleading colour is 



(5.29) 



^0 



+ 



-4+- + 2Z 
z 
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1 r 
1 

+- 
e 



-12 - 10H(0; z) - 8H(1; z) + - (16 + 8H(0; z) + 8H(1; z)) +z(^ + 5H(0; z) + 4H(1; z)) 
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-— + 47r^ - 22H(0; z) - 24H(1; z) - 22H(0, 0; z) - 20H(0, 1; z) - 16H(1, 0; z) 



+ 



- 16H(1, 1; z) + - (67 - + 32H(0; z) + 32H(1; z) + 16H(0, 0; z) + 16H(0, 1; z) 

Z o 

+ 16H(1, 0; z) + 16H(1, 1; z)) + z(y - It? + yH(0; z) + 7H(1; z) + 11H(0, 0; z) 
+ 10H(0, 1; z) + 8H(1, 0; z) + 8H(1, 1; 2)) 

-184 + ^ + 1^ - ^TT^HC-l; z) - 1^H(0; + %t?^{Q, z) - 91H(1; 

+ 87r^H(l; z) + 8H(-1, 0; z) - 38H(0, 0; z) - 44H(0, 1; z) - 44H(1, 0; z) 

- 48H(1, 1; z) - 16H(-1, -1, 0; z) + 8H(-1, 0, 0; z) + 16H(0, -1, 0; z) - 46H(0, 0, 0; z) 

- 44H(0, 0, 1; z) - 40H(0, 1, 0; z) - 40H(0, 1, 1; z) - 16H(1, 0, 0; z) - 32H(1, 0, 1; z) 



-24H(l,l,0;z)-32H(l,l,l;z) + -( 

Z Zi 



1.485 407r2 8OC3 4 



7r2H(-l;z) + 134H(0;z) 



20 

- — 7r^H(0; z) + 134H(1; z) - 87r^H(l; z) + 64H(0, 0; ^) + 64H(0, 1; z) + 64H(1, 0; z) 

+ 64H(1, 1; z) - 16H(-1, -1, 0; z) + 8H(-1, 0, 0; z) + 32H(0, 0, 0; z) + 32H(0, 0, 1; z) 
+ 32H(0, 1, 0; z) + 32H(0, 1, 1; z) + 16H(1, 0, 0; z) + 32H(1, 0, 1; z) + 24H(1, 1, 0; z) 
97 It? 



+ 32H(l,l,l;z)) 



^ - ^7r2H(-l; z) + 11H(0; z) - \t?^{^; z) 

+ 11H(1; 2) - 47r2H(l; ^) + 8H(-1, 0; z) + yH(0, 0; z) + 19H(0, 1; z) + 16H(1, 0; z) 

+ 14H(1, 1; z) - 8H(-1, -1, 0; z) + 4H(-1, 0, 0; z) + 23H(0, 0, 0; z) + 22H(0, 0, 1; z) 
+ 20H(0, 1, 0; z) + 20H(0, 1, 1; z) + 8H(1, 0, 0; z) + 16H(1, 0, 1; z) + 12H(1, 1, 0; z) 

+ 16H(l,l,l;z)) 



+ 0(e) 



(5.30) 



The four-quark antenna function contains two different quark flavours. Since it is 
symmetric under the interchange of the quark-antiquark pairs, we have to consider only- 
one crossing, resulting in {q,qq'q'), which integrates to 



'^q,qq'<l' 



+ - 
+ - 

+ 
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+ ^H(0;z) 


e 





50 

T 
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"9" 



1^ 



-3H(l;z)) 

505 IOtt^ 100,,, , 100,,, , 16,,, , 16,,, , 16„, 

^ - + -^H(0; z) + —nil; z) + y H(0, 0; z) + -H(0, l;z) + -H(l, 0; z) 



16,,,^ -, X 1 / 505 IOtt^ 
+ -H(l,l;.) + -(-- + ^ 



100 



H(0;z) 



100 



H(l;z) 



16 



H(0,0;z) 
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- y H(0, 1; z) - -H(l, 0; z) - y H(l, 1; z)) + ^(--- + _ 

- ^H(0; z) - ^H(l; - ^H(0, 0; z) - ^H(0, 1; z) - ^H(l, 0; z) - ^H(l, 1; z)) 



+ 0{e) 



(5.31) 



The quark-initiated one-loop virtual gluon-gluon antenna of the form (g, qg) have a 
leading and subleading colour contribution and a quark loop piece. Their integrals are: 



+- 



+- 



1 / 3 



3 1, 



+ H(0; z) + -H(l; z) + -{--- H(0; z) - -H(l; z)) -z{- + -H(0; z) + -H(l; z)) 



8 2 



53 TT 2 5 3 

■- - — + -H(0; + -H(l; z) + 4H(0, 0; z) + 2H(0, 1; z) + 2H(1, 0; z) + -H(l, 1; z) 
3d 3 3 4 2 

1 79 TT^ 13 

+ - + y - -g-H(0; - 2H(1; z) - 4H(0, 0; - 2H(0, 1; z) - 2H(1, 0; z) 
- ^H(l, 1; z)) + + ^ + iH(0; - ^H(l; z) - 2H(0, 0; - H(0, 1; z) 



H(l,0;z)--H(l,l;2)) 



+ 



r 571 357r2 



5C3 - ^H(0; z) - ^7r2H(0; z) + 4h(1; ^) - 7^'H(1; z) - H(0, 0; z) 

3d 



Z5i 

36 



11 



108 12 

+ ^H(0, 1; z) - ^H(l, 0; z) + ^H(l, 1; 2) + 10H(0, 0, 0; z) + 6H(0, 0, 1; z) + 4H(0, 1, 0; z) 
+ 4H(0, 1, 1; z) + 5H(1, 0, 0; z) + 4H(1, 0, 1; z) + 4H(1, 1, 0; z) + ^H(l, 1, 1; z) 



1 ,3733 377r2 43,,, , 11 n.w x 25,,, , 3 n.w x 

+ 7 + -PT + - ^H(0; z) + -7r2H(0; ^) - -H(l; z) + -tt^HII; z) 



z M32 12 

14 19 9 

- 6H(0, 0; z) - yH(0, 1; z) - -^H(l, 0; z) - -H(l, 1; z) - 10H(0, 0, 0; z) 

- 6H(0, 0, 1; z) - 4H(0, 1, 0; z) - 4H(0, 1, 1; z) - 5H(1, 0, 0; z) - 4H(1, 0, 1; z) 

- 4H(1, 1, 0; z) - ^H(l, 1, 1; ^)) + + ^ + ^ + |h(0; z) + ^tt^HCO; z) 
+ y H(l; z) + ^7r2H(l; z) + ^H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z) - ^H(l, 1; z) 

- 5H(0, 0, 0; z) - 3H(0, 0, 1; 2) - 2H(0, 1, 0; z) - 2H(0, 1, 1; z) - ^H(l, 0, 0; z) 



- 2H(1, 0, 1; z) - 2H(1, 1, 0; z) - -H(l, 1, 1; z)) 

+ 0{e) , 



(5.32) 



q1 _ 



1 ri 

5 



1 

2i 



+ 



L4 



+ H(0;z) + -H(l;z)--(2 + H(0;z) 



+ 2^(1;^)) 



z(UiH(0;z) 



+ 4H(1;^))_ 
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+ 



1 



+ ^ (-^ + ^ - 4H(0; z) - 2H(1; - 2H(0, 0; z) - H(0, 1; z) - H(l, 0; z) 
- ^H(l, 1; z)) + z{-l + ^ - ^H(0; z) - ^H(l; z) - H(0, 0; z) - 1h(0, 1; z) 



+ 



51 

T 



iH(l,0;^)-^H(l,l;z)) 



+ ^H(0, 1; z) + 2H(1, 0; + ^H(l, 1; + 4H(0, 0, 0; z) + 2H(0, 0, 1; z) + 2H(0, 1, 0; z) 

+ H(0, 1, 1; z) + H(l, 0, 0; z) + H(l, 0, 1; z) + ^H(l, 1, 1; 2;) + i (-1^ + tt^ + 1^ 

31 1 67 7 

- -H(0; z) + -7r2H(0; z) - — H(l; z) + y^vt^HCI; z) - 8H(0, 0; z) - 4H(0, 1; z) 

- 4H(1, 0; z) - 2H(1, 1; - 4H(0, 0, 0; z) - 2H(0, 0, 1; z) - 2H(0, 1, 0; z) - H(0, 1, 1; 

- H(l, 0, 0; z) - H(l, 0, 1; z) - ^H(l, 1, 1; z)) + z{-2 + ^ + ^ - 2H(0; 



1 



TT^BiO-z) -R{l;z) 



\^R{l;z)-'^-R{0,0;z) 



24 



-H(0,l;z) -H(l,0;z) 



4 

^H(l, 1; z) - 2H(0, 0, 0; z) - H(0, 0, 1; z) - H(0, 1, 0; z) - ^H(0, 1, 1; z) 



1 



1 



1, 



-H(l, 0, 0; z) - -H(l, 0, 1; z) - -H(l, 1, 1; z)) 



+ 0(e), 



(5.33) 



qI 



1 

+- 
e 



TT^ 70 2Q 1 1 14 

" T " 36^^°' " 36^^^' ~ °' " 3"^°' ^' " 3^^^' °' + 27 

- Y - ^H(0; z) - ^H(l; - ^H(0, 0; z) - l;z)- ^H(l, 0; z)) 

+ 0{e) . (5.34) 

6. Gluon initiated antennae 

The gluon-initiated initial-final antenna functions are obtained from final-final antenna 
functions listed in [24] by crossing a gluon into the initial state. As in the quark case, 
their unintegrated forms are pure kinematical crossings of the final- final expressions, with 
no symmetry factors or polarization sums multiplied on them. Also, we do not provide a 
decomposition into sub-antennae here. The precise decomposition depends on the require- 
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ments and symmetries of the process under consideration, and is normally performed in 
the context of an actual calculation. 

As in the quark-initiated case, some attention has to be paid in the crossing of the 
virtual one-loop antenna functions from the final-final kinematics [24] to the initial-final 
kinematics. This crossing requires the analytical continuation of the poly logarithmic func- 
tions to the relevant phase space, again consisting of three different Riemann sheets [53,54]. 

All gluon- initiated initial- final antenna functions are listed in Table 2, their integrated 
forms are collected in the following. 



6.1 Quark-quark antennae 

The NLO quark-antiquark antenna function to order is: 

= 



1 

+- 

e 



+ 



1 - 2z + 2r 
(l-2z + 2z2)(H(0;z) + H(l;z)) 



2 + (1 - 2z + 2z2) + H(0, 0; z) + H(0, 1; z) + H(l, 0; z) + H(l, 1; z)) 



4 + 2H(0; z) + 2H(1; z) + {l - 2z + 2z^) (-^ - ^7r2H(0; z) - ^7r2H(l; z) 

+ H(0, 0, 0; z) + H(0, 0, 1; z) + H(0, 1, 0; z) + H(0, 1, 1; z) + H(l, 0, 0; z) + H(l, 0, 1; z) 

+ H(l,l,0;z)+H(l,l,l;z))" 

+ 0(e3). (6.1) 

For the NNLO double real radiation antenna functions, we have the leading and subleading 
colour contributions to {g,gqq), whose integrals read 

1 r 



1 



1 



3 - 62; Qz^ 



+ 



L 4 



13 4 7 

— + -H(0; z) + 6H(1; z) + z(-15 - 19H(0; z) - 12H(1; z)) 



A9 



+- 



+ ^2 (_ + 12H(0; z) + 12H(1; z)) 

173 
12 



27 



13 



11, 



- 27r^ + — H(0; z) + — H(l; z) - 2H(-1, 0; z) + — H(0, 0; z) + 7H(0, 1; z) 



+ 



1 8 8 8 9 227r2 

+ 11H(1, 0; z) + 12H(1, 1; + - (- - -H(l; z))+z{-- + ^- 16H(0; z) 

- 30H(1; z) - 4H(-1, 0; z) - 43H(0, 0; z) - 38H(0, 1; z) - 22H(1, 0; z) - 24H(1, 1; z)) 
+ ( Y - ^ + ^H(0; z) + y H(l; z) - 4H(-1, 0; z) + 26H(0, 0; z) + 24H(0, 1; z) 
+ 22H(1,0;2) +24H(l,l;z)) 

- ^ - 41C3 + 7r2H(-l; .) + ^H(0; .) - ^Tr'm z) + ^R{1; z) - 57r'R{l; z) 
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+ 2H(-1, 0; z) + — H(0, 0; z) + yH(0, 1; z) + 13H(1, 0; z) + 13H(1, 1; z) 

19 

+ 4H(-1, -1, 0; z) - 8H(-1, 0, 0; z) - 4H(-1, 0, 1; z) - 8H(0, -1, 0; z) + yH(0, 0, 0; z) 

+ 11H(0, 0, 1; z) + 10H(0, 1, 0; z) + 14H(0, 1, 1; z) + 17H(1, 0, 0; z) + 22H(1, 0, 1; z) 

2 , , 1 , 502 27r2 2 , , 

+ 22H(1, 1, 0; z) + 24H(1, 1, 1; z) + ^(-1, ^) + - ("^ + -3- - gHlO; ^) 

+ ^H(l; .) + ^H(-l, 0; .) - ^H(l, 0; .) - ^H(l, 1; .)) + .(-iH^ + 1^ + 86(3 

917 191 8Q 

+ 27r2H(-l; z) - =^H(0; z) + ^7r2H(0; z) - ^H(l; z) + 107r2H(l; z) 

56 116 

+ — H(-l, 0; z) H(0, 0; z) - 32H(0, 1; z) - 62H(1, 0; z) - 60H(1, 1; z) 

3 3 

+ 8H(-1, -1, 0; z) - 16H(-1, 0, 0; z) - 8H(-1, 0, 1; z) - 91H(0, 0, 0; z) - 86H(0, 0, 1; z) 

- 68H(0, 1, 0; z) - 76H(0, 1, 1; z) - 34H(1, 0, 0; z) - 44H(1, 0, 1; z) - 44H(1, 1, 0; z) 

- 48H(1, 1, 1; z)) + _ ^ _ 92^3 + 27r2H(-l; z) - ^H(0; z) - ^7r2H(0; ^) 

1 89 44 398 984 

+ — H(l; z) - 107r2H(l; z) + -H(-l, 0; z) + — H(0, 0; z) + — H(0, 1; z) 

196 196 

+ — H(l,0;z) + — H(l,l;z) +8H(-l,-l,0;z) -16H(-l,0,0;z) -8H(-l,0,l;z) 

- 16H(0, -1, 0; z) + 54H(0, 0, 0; z) + 52H(0, 0, 1; z) + 40H(0, 1, 0; z) + 48H(0, 1, 1; z) 
+ 34H(1, 0, 0; z) + 44H(1, 0, 1; z) + 44H(1, 1, 0; z) + 48H(1, 1, 1; z)) 

+ 0(e) (6.2) 
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+- 
e 



4 - 8z + 8z- 
13 

T 

+ ^2(12 + 16H(0; z) + 16H(1; z)) 



+ 7H(0; z) + 8H(1; z) - z(l4 + 14H(0; z) + 16H(1; z)) 



+ 



39 Ttt^ 15 

- - — + — H(0; z) + 13H(1; z) + 15H(0, 0; z) + 14H(0, 1; z) + 18H(1, 0; z) 

Zi d ^ 

75 147r^ 

+ 16H(1, 1; z) + + ^ 32H(0; z) - 28H(1; z) - 30H(0, 0; z) - 28H(0, 1; z) 

- 36H(1, 0; z) - 32H(1, 1; z)) + (38 - + 24H(0; z) + 24H(1; z) + 36H(0, 0; z) 
+ 32H(0, 1; z) + 36H(1, 0; + 32H(1, 1; z)) 

301 _ 9^ _ 8|3 ^ ^,.„(_1, + i??H(0; .) - f .^H(0; .) + 39H(1; .) - f .^H(l; .) 

27 

- 24H(-1, 0; z) + — H(0, 0; z) + 15H(0, 1; z) + 30H(1, 0; z) + 26H(1, 1; z) 
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+ 16H(-1, -1,0; z) - 8H(-1, 0, 0; z) - 16H(0, -1, 0; z) + 31H(0, 0, 0; z) + 30H(0, 0, 1; z) 
+ 28H(0, 1, 0; z) + 28H(0, 1, 1; z) + 26H(1, 0, 0; z) + 36H(1, 0, 1; z) + 28H(1, 1, 0; z) 

+ 32H(1, 1,1; z) + ^H(-l, 0; z) + - (- - -H(0; z)) + ^(-99 + + ^ 

+ ^7r2H(-l; z) - ^H(0; + ^n^m ^) " 75H(1; + y tt^HCI; - ^H(-l, 0; z) 

- — H(0, 0; z) - 64H(0, 1; z) - 60H(1, 0; z) - 56H(1, 1; z) + 32H(-1, -1, 0; z) 

- 16H(-1, 0, 0; z) - 62H(0, 0, 0; z) - 60H(0, 0, 1; z) - 56H(0, 1, 0; z) - 56H(0, 1, 1; z) 

- 52H(1, 0, 0; z) - 72H(1, 0, 1; z) - 56H(1, 1, 0; z) - 64H(1, 1, l\z)) + z^ (98 - IOtt^ 
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^ + -7r"H(-l; + 76H(0; z) - y7r^H(0; z) + 76H(1; z) - y7r2H(l; 

+ 48H(0, 0; z) + 48H(0, 1; z) + 48H(1, 0; z) + 48H(1, 1; z) + 16H(-1, -1, 0; z) 
- 8H(-1, 0, 0; z) - 16H(0, -1, 0; z) + 76H(0, 0, 0; z) + 72H(0, 0, 1; z) + 64H(0, 1, 0; z) 

+ 64H(0, 1, 1; z) + 68H(1, 0, 0; z) + 72H(1, 0, 1; z) + 64H(1, 1, 0; z) + 64H(1, 1, 1; z)) 

+ 0(e). (6.3) 



The one-loop corrections to the (g, qq) antenna function have a leading and subleading 
colour term as well as a closed quark loop contribution. These integrate to 
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2 

^ ^-H(0; z) - hl{l; z) + z{^ + 3H(0; z) + H(l; z)) 
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H(l; z) - 8H(0, 0; z) - 3H(0, 1; z) - 2H(1, 0; z) - H(l, 1; z)) 



+ 



+ z2(---3H(0;z)-H(l;z)) 
-4 + y - y H(0; z) - ^H(l; z) - 4H(0, 0; - ^H(0, 1; z) - H(l, 0; z) - ^H(l, 1; z) 
+ - TT^ + ^H(0; z) + ^H(l; z) + 8H(0, 0; z) + 3H(0, 1; z) + 2H(1, 0; z) + H(l, 1; z)) 
+ .^(7r^-HH(0;.)-H 

- y H(0, 1; z) - y H(l, 0; z) - ^H(l, 1; z) - 9H(0, 0, 0; z) - 4H(0, 0, 1; z) - 2H(0, 1, 0; z) 

- -H(0, 1, 1; z) - H(l, 0, 0; z) - H(l, 0, 1; z) - -H(l, 1, 1; 2) + z(- - ^ - ^ 

+ Jh(0; - ^7r2H(0; z) + iH(l; z) - lir^m- ^) + ^H(0, 0; z) + ^H(0, 1; z) 
ZD 4 6 3 3 

14 11 

+ — H(l, 0; z) + — H(l, 1; z) + 18H(0, 0, 0; z) + 8H(0, 0, 1; z) + 4H(0, 1, 0; z) 
o o 

+ 3H(0, 1, 1; z) + 2H(1, 0, 0; z) + 2H(1, 0, 1; z) + H(l, 1, 1; z)) + + ^ 
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+ ^7r2H(0; z) + ^7r2H(l; z) - y H(0, 0; z) - y H(0, 1; z) - ^H(l, 0; z) 

- ^H(l, 1; z) - 18H(0, 0, 0; z) - 8H(0, 0, 1; z) - 4H(0, 1, 0; z) - 3H(0, 1, 1; z) 

- 2H(1, 0, 0; z) - 2H(1, 0, 1; z) - H(l, 1, 1; z)) 
+ 0{e) 

and 



(6.4) 



1 



-^H(0; z) - ^H(l; ^) + z(H(0; ^) + H(l; z)) - z2(H(0; z) + H(l; z)) 



Ir7r2 



+- 
e 



■— - ^H(0; ^) - ^H(l; z) - 2H(0, 0; z) - ^H(0, 1; z) - 2H(1, 0; z) - ^H(l, 1; z) 



+ ^(~7 ~ ^ + 4H(0, 0; z) + 3H(0, 1; z) + 4H(1, 0; z) + 3H(1, 1; z)) 



+ z'^(—- 4H(0, 0; z) - 3H(0, 1; z) - 4H(1, 0; z) - 3H(1, 1; z)) 
6 



+ 



5 TT 

L4 ~ 12 

3 



^ _ 3H(0; z) + ^7r2H(0; z) - 4H(1; z) + ^7r2H(l; - ^ 



12 



H(0,0;z) 



2H(0, 1; z) - ^H(l, 0; z) - ^H(l, 1; z) - 5H(0, 0, 0; z) - 4H(0, 0, 1; z) 

- 4H(0, 1, 0; z) - ^H(0, 1, 1; z) - 4H(1, 0, 0; z) - 4H(1, 0, 1; z) - 3H(1, 1, 0; z) 

- Ih(1, 1, l-z)+z{-^ + ^ + 3C3 - ^H(0; z) - ^n^m z) - ^H(l; z) - y ^'H(1; z) 

+ H(0, 0; z) + H(l, 0; z) + 10H(0, 0, 0; z) + 8H(0, 0, 1; z) + 8H(0, 1, 0; z) + 7H(0, 1, 1; z) 
+ 8H(1, 0, 0; z) + 8H(1, 0, 1; z) + 6H(1, 1, 0; z) + 7H(1, l,hz))+z^ (-(s + H(0; z) 

+ ^7r2H(0; z) + H(l; z) + ^7r2H(l; z) - 10H(0, 0, 0; z) - 8H(0, 0, 1; z) - 8H(0, 1, 0; z) 



- 7H(0, 1, 1; z) - 10H(1, 0, 0; z) - 8H(1, 0, 1; z) - 8H(1, 1, 0; z) - 7H(1, 1, 1; z)) 
+ 0{e) 
and 



(6.5) 



+- 



-1 2z 2^2 
-3 ~ Y ~3~ 

iH(0; z) + ^H(l; z) + z(-^H(0; z) - ^H(l; z)) + z^ (^H(0; z) + ^H(l; z)) 



2 2, 



^ - ^ + ^H(0, 0; z) + iH(0, 1; z) + ^H(l, 0; z) + ^H(l, l-z)+z{—~ |H(0, 0; z) 
- ^H(0, 1; z) - ^H(l, 0; z) - ^H(l, 1; z)) +z'{-^ + ^H(0, 0; 2) + ^H(0, 1; 2) 
+ ^H(l,0;z) + ^H(l,l;z)) 
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+ 0{e) 



(6.6) 



6.2 Quark-gluon antennae 

The integrated gluon-initiated {g,qq) antenna function at NLO is through to order e^: 



V" = 

9,91 



1 

+ - 
e 



+ 



\ + '^6{l -z)- Vo{z) ---2z + 2z^ 



-\ + '^(1 - ^)(^ - x) - i^o(^) + ^i(^) + ^"(^5 ^) + ^H(l; z) - Y^H(0; z) 
_ i + H(0; z) + H(l; z)) + {-z + z^) (2H(0; z) + 2H(1; z)) 

1 - — + <5(1 - z)(- - — - ^) + (-- + -)Po(^) + - o^2(^) - :^H(0;z) 



4 

iH(l; z) + |h(0, 0; z) + |h(0, 1; z) + |h(1, 0; z) + |h(1, 1; z) 



^ :(7H(0;.) 



1 -zM 



+ H(0, 0; z) + H(0, 1; z) + H(l, 0; z)) + ^ (- J + ^ - ^H(0; z) - ^H(l; z) 

J2 



- H(0, 0; z) - H(0, 1; z) - H(l, 0; z) - H(l, 1; z)) + (z - z^) (y - 2H(0, 0; z) 

- 2H(0, 1; z) - 2H(1, 0; z) - 2H(1, 1; z)) 



vr 



3 



35C3 



77r2 TT^ 



7C3^ 



37r^ 



7C3> 



7 ^ 



-^)Pi(. 



I?2(^) + \v^{z) + H(0; z) - ^7r2H(0; z) + H(l; z) - ^7r2H(l; z) - ^H(0, 0; z) 



- iH(0, 1; z) - ^H(l, 0; z) - ^H(l, 1; z) + ^H(0, 0, 0; z) + ^H(0, 0, 1; z) + ^H(0, 1, 0; z) 
+ ^H(0, 1, 1; z) + ^H(l, 0, 0; z) + ^H(l, 0, 1; z) + ^H(l, 1, 0; z) + ^H(l, 1, 1; z) 



+ 



1 



I- z 



■^H(0; z) + ^7r2H(0; z) - ^H(0, 0; z) - ^H(0, 1; z) - ^H(l, 0; z) 



H(0, 0, 0; z) - H(0, 0, 1; z) - H(0, 1, 0; z) - H(0, 1, 1; z) - H(l, 0, 0; z) - H(l, 0, 1; z) 



7C3 7^ 



- H(l, 1, 0; z)) +-(-- + — + ^ - -H(0; z) + -7r2H(0; z) - -H(l; z) + -7r^R{l; z) 

- ^H(0, 0; z) - ^H(0, 1; z) - ^H(l, 0; z) - ^H(l, 1; z) - H(0, 0, 0; z) - H(0, 0, 1; z) 

- H(0, 1,0; z) - H(0, 1, 1; z) - H(l, 0, 0; z) - H(l, 0, 1; z) - H(l, 1, 0; z) - H(l, 1, 1; z)) 

+ - '^') + ^^'H(0; z) + j7r2H(l; z) - 2H(0, 0, 0; z) - 2H(0, 0, 1; z) - 2H(0, 1, 0; z) 



- 2H(0, 1, 1; z) - 2H(1, 0, 0; z) - 2H(1, 0, 1; z) - 2H(1, 1, 0; z) - 2H(1, 1, 1; z)) 
+ 0(e3) . 



(6.7) 
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In crossing one of the gluons in the final-final quark-gluon antenna function into 
the initial state, one has to distinguish two cases, depending on whether the gluon is colour- 
connected to the quark or not. If the crossed gluon is colour-connected to the quark, we 
obtain the following integrated antenna function: 



1 



1 

+3 



36(1 - z) 



1 O f\ 

16 + —(5(1 -z)- 6Vo{z) Uz + \^z^ 

3 z 



+ 



1 



85 



1429 ItP- 



38 



e2 L 4 



+ 5(1 - z) {— -)- -Vo{z) + l2V,{z) + 26H(0; z) + 32H(1; z) - — 



12H(0; z) 



+ z{-^ - 34H(0; z) - 28H(1; z)) + ^(-y - 12H(0; z) - 12H(1; z)) 
+ ^2(31 + 28H(0; z) + 28H(1; z)) 

223 2l7r2 ,,23959 18l7r2 , , 1429 Mtt^, , , 76 , , 

+ - '^^^2 ~ - + + -^)^o(^) + ^V,{z) - 12V2{z) 



36 



36 



4 2 
595 85 

+ -n7H(0; z) + — H(l; z) - 2H(-1, 0; z) + 49H(0, 0; z) + 52H(0, 1; z) -h 65H(1, 0; z) 



12 



1 



,76 



+ 64H(1, 1; z) - (^H(0; z) + 24H(0, 0; z) + 24H(0, 1; z) + 24H(1, 0; z)) 

-L ,2 o 

+ ^ (-^ + ^ - f H(0; z) - ^H(l; z) - 24H(0, 0; z) - 24H(0, 1; z) - 24H(1, 0; z) 

- 24H(1, 1; z))+z(-'^ + ^ - — H(0; 2) - 61H(1; z) - 4H(-1, 0; z) - 76H(0, 0; z) 

' ^ 4 3 6 

- 68H(0, 1; z) - 58H(1, 0; z) - 56H(1, 1; z)) + 



;^-llvr2 + ^H(0;z) + 62H(l;.) 



4H(-1, 0; z) + 62H(0, 0; + 56H(0, 1; z) 58H(1, 0; z) + 56H(1, 1; z)) 



+ 



7793 13857r2 479(3 ,,389623 68577r2 1637r4 1027(3^ 
^ +5(1 - z)( \ —) 



18 



36 72 3 ^ 'V 2592 432 

, 23959 18l7r2 x_ , , ,1429 287r2. , , 76_ , , . , „ , , 

+ ("^16" + + 50C3)25o(^) + {—^ —)-Diiz) - -V2iz) + 8^3(2;) + 37r2H(-l; ^) 



18 



18 



2801 



223 



317, 



+ H(0; z) - 247r^H(0; z) + — H(l; z) - 277r^H(l; z) + 10H(-1, 0; z) + — H(0, 0; z) 
60 2 4 

595 

-^H(0, 1; z) + 78H(1, 0; z) + 85H(1, 1; z) + 28H(-1, -1, 0; z) - 20H(-1, 0, 0; z) 

- 4H(-1, 0, 1; z) - 44H(0, -1, 0; z) + 95H(0, 0, 0; z) + 98H(0, 0, 1; z) + 94H(0, 1, 0; z) 
+ 104H(0, 1, 1; z) + 101H(1, 0, 0; z) + 130H(1, 0, 1; z) + 118H(1, 1, 0; z) + 128H(1, 1, 1; z) 



1 ,5^2 
+ z -^ + 16C3 



i^H(0; .) + ^n^m z) - ipH(0, 0; .) - ipH(0, 1; .) 



l-z^ 6 
137 

— H(l, 0; z) + 8H(0, -1, 0; z) - 48H(0, 0, 0; z) - 48H(0, 0, 1; z) - 44H(0, 1, 0; z) 



-48H(0,l,l;z) -44H(l,0,0;z) - 48H(1, 0, 1; z) - 48H(1, 1, 0; z)) + z{- 



9395 
36 
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70l7r2 376(3 , 2tt/ . x 6349 
+ + + 47r2H(-l; - ^ 



H(0; z) + ^7r2H(0; 2) - ^H(l; z) 



265, 

T 



36 3 

201 307 
-2xjn . .A .jj(o, 0; z) - — H(0, 1; z) - 130H(1, 0; z) - 122H(1, 1; z 



+ 247r2H(l;z) ^ 3 

+ 32H(-1, -1, 0; z) - 28H(-1, 0, 0; z) - 8H(-1, 0, 1; z) - 160H(0, 0, 0; z) 

- 152H(0, 0, 1; z) - 128H(0, 1, 0; z) - 136H(0, 1, 1; z) - 94H(1, 0, 0; z) - 116H(1, 0, 1; 2) 

\, 46025 15l7r2 2 , , 

- 104H(1, 1, 0; z) - 112H(1, 1, 1; z)) + - ( — + + 54(3 + o7r^H(-l; z) 



216 



9 



^^^^H(0;z) + ^7r2H(0;z)- 



1 823 68 1 Q6 

-^H(l; z) + IOtt^HCI; z) + -H(-l, 0; z) - — H(0, 0; z) 



18 ' ' ' 3 

196 239 248 

- — H(0, 1; z) - — H(l, 0; z) - — H(l, 1; z) + 8H(-1, -1, 0; z) - 4H(-1, 0, 0; z) 

+ 8H(0, -1, 0; z) - 48H(0, 0, 0; z) - 48H(0, 0, 1; z) - 44H(0, 1, 0; z) - 48H(0, 1, 1; z) 

- 36H(1, 0, 0; z) - 48H(1, 0, 1; z) - 44H(1, 1, 0; z) - 48H(1, 1, 1; z)) + z^ (-*^^^^ ^^^^ 



27 



18 



484(3 



10 



862 



1130, 



70 



+ — 7r^H(-l; z) + -^H(0; z) - 287r^H(0; z) + -r-H(l; z) - — 7r^H(l; z) 



44 460 

+ — H(-l, 0; z) + 168H(0, 0; z) + H(0, 1; z) + 124H(1, 0; z) + 124H(1, 1; z) 

3 3 

+ 24H(-1, -1, 0; z) ~ 24H(-1, 0, 0; z) - 8H(-1, 0, 1; z) - 32H(0, -1, 0; z) 
+ 130H(0, 0, 0; z) + 124H(0, 0, 1; z) + 104H(0, 1, 0; z) + 112H(0, 1, 1; z) 

+ 102H(1, 0, 0; z) + 116H(1, 0, 1; z) + 108H(1, 1, 0; z) + 112H(1, 1, 1; z)) 

+ Oie), 



(6.8) 



while the crossing the middle gluon (which is colour-connected to the other two gluons) 
results in: 



T,0 ^ 

9', 991 



1 r 
1 



4^(1 



18 + 6(5(1 -z)- 8Vo{z) ---12z + 12z^ 



1 + ,(1 _ ,)(Z^ _ 1^) _ i2Po(.) + 16P,(.) + 26H(0;.) - '-^^ + i(-28 



2 ' 'M 

- 16H(0; z) - 16H(1; z)) + 36H(1; z) + z(-29 - 40H(0; z) - 24H(1; z)) 
+ z2(34 + 24H(0; z) + 24H(1; z)) 



1 

+- 
e 



116(3, / 75 207r2 
^) + ( \ 



)X'o(2) + 24Pi(z)-16I>2(-2) 



85 407r2 ^.417 ^ 2 

____ + ,(l_,)(__5.^ . . . 

133 

+ ^H(0; z) + 45H(1; z) - 16H(-1, 0; z) + 56H(0, 0; z) + 52H(0, 1; z) + 70H(1, 0; z) 
+ 72H(1, 1; + (-24H(0; z) - 36H(0, 0; z) - 32H(0, 1; - 32H(1, 0; z)) 



+ 



1 



,27r^ 



1 . 785 



116, 



^ ^ ^ (— + 8H(-1, 0; z) - 4H(0, 0; z)) +-{-— + 61:^ - — H(0; - 56H(1; z) 
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59 

- 8H(-1, 0; z) - 32H(0, 0; z) - 32H(0, 1; z) - 32H(1, 0; z) - 32H(1, 1; z)) + z{- — 
467r^ 

+ — 29H(0; z) - 58H(1; z) - 8H(-1, 0; z) - 92H(0, 0; z) - 80H(0, 1; z) 

- 44H(1, 0; z) - 48H(1, 1; z)) + ^^(^ - ^ + — H(0; z) + 68H(1; z) 

- 8H(-1, 0; z) + 52H(0, 0; z) + 48H(0, 1; z) + 44H(1, 0; z) + 48H(1, 1; z)) 
3493 112l7r2 .,2275 1257r2 GItt^ 417 



+ - - + - ^' - — + -90- - «6) + "°(^) (-^ + 

+ + (75 - ^)Vi{z) - 2iV,{z) + + f ^'^H(-li^) + HIh(0;j) 

89 8fi 1055 

- — 7r2H(0; z) + 85H(1; z) - z) + -7-H(0, 0; z) + 133H(0, 1; z) 

o O D 

262 

+ H(l, 0; z) + 90H(1, 1; z) + 40H(-1, -1, 0; z) - 68H(-1, 0, 0; z) - 32H(-1, 0, 1; z) 

3 

- 36H(0, -1, 0; z) + 116H(0, 0, 0; z) + 112H(0, 0, 1; z) + 84H(0, 1, 0; z) + 104H(0, 1, 1; z) 
+ 146H(1, 0, 0; z) + 140H(1, 0, 1; z) + 148H(1, 1, 0; z) + 144H(1, 1, 1; z) 

1 52 166 

+ (+24C3 - 79H(0; z) + — 7r2H(0; z) - — H(0, 0; z) - 48H(0, 1; z) - 48H(1, 0; z) 

- 76H(0, 0, 0; z) - 72H(0, 0, 1; z) - 56H(0, 1, 0; z) - 64H(0, 1, 1; z) - 64H(1, 0, 0; z) 

- 64H(1, 0, 1; z) - 64H(1, 1, 0; z)) + ^ (+28(3 - y tt^HC-I; z) + ^7r2H(0; z) 

- 32H(-1, -1, 0; + 40H(-1, 0, 0; z) + 16H(-1, 0, 1; z) + 24H(0, -1, 0; z) 

- 12H(0, 0, 0; z) - 8H(0, 0, 1; z)) + i (-^ + ^ + y 7r2H(-l; - ^-^m z) 

785 38 20 232 

+ 167r2H(0; z) - — H(l; z) + yTT^HCl; z) - yH(-1, 0; - — H(0, 0; z) 

232 314 

- — H(0, 1; z) - — H(l, 0; z) - 112H(1, 1; z) + 24H(-1, -1, 0; z) - 36H(-1, 0, 0; z) 

- 16H(-1, 0, 1; 2) - 24H(0, -1, 0; z) - 64H(0, 0, 0; z) - 64H(0, 0, 1; z) - 56H(0, 1, 0; z) 

- 64H(0, 1, 1; z) - 72H(1, 0, 0; z) - 64H(1, 0, 1; z) - 68H(1, 1, 0; z) - 64H(1, 1, 1; z) 

256C3X / 3571 3777r2 , 10 , 857^,, , 128 n,,, , 

+ -3^) + + + ^^^^3 + Y^'H(-1; z) - — H(0; z) + — 7r2H(0; z) 

- 59H(1; z) + ^7r2H(l; z) + 36H(-1, 0; z) - ^H(0, 0; z) - 58H(0, 1; z) 
344 

- — H(l, 0; z) - 116H(1, 1; z) + 8H(-1, -1, 0; z) - 28H(-1, 0, 0; z) - 16H(-1, 0, 1; z) 

- 196H(0, 0, 0; z) - 184H(0, 0, 1; z) - 144H(0, 1, 0; z) - 160H(0, 1, 1; z) - 76H(1, 0, 0; z) 

- 88H(1, 0, 1; z) - 92H(1, 1, 0; z) - 96H(1, 1,1; z))+z^ {— 184(3 

8 80 560 88 

+ 47r2H(-l; z) + -H(0; z) - -ir^BiO; z) + — H(l; z) - 207r2H(l; z) + yH(-1, 0; z) 

584 

+ 224H(0, 0; z) + — H(0, 1; z) + 136H(1, 0; z) + 136H(1, 1; z) + 16H(-1, -1, 0; z) 



2 
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- 32H(-1, 0, 0; z) - 16H(-1, 0, 1; z) ~ 32H(0, -1, 0; z) + 108H(0, 0, 0; z) 
+ 104H(0, 0, 1; z) + 80H(0, 1, 0; z) + 96H(0, 1, 1; z) + 68H(1, 0, 0; z) + 88H(1, 0, 1; z) 

+ 88H(1, 1, 0; z) + 96H(1, 1, 1; z)) 

+ 0{e). (6.9) 



Crossing the gluon in the E'^-type antenna functions into the initial state results in 
{g, qq' q')-antemiSL at leading and subleading colour, which integrate to: 



cO 



1 



1 r 3 19 



- o - - ^) + i^o(^) + H(0; ^) + ^ + ^( 9 + H(0; ^)) - ^ 



+- 



2 18 

7 7r2 



373 Stt^ 19_ , , 4_ , , 5„,^ , 4 
3 3 +^(l-^H-io8+l^) + ¥^°(^)-3^^(^)-6^(«^^) + 3(r3i) 



H(0;2:) 



+ 



1 52 8 8 13 71"^ 

- 3H(1; z) + 3H(0, 0; z) + 2H(0, 1; ^) + - (- + -H(0; z) + -H(l; ^)) + z(y - y 

25 16 16 

+ — H(0; z) + 3H(1; z) + 3H(0, 0; z) + 2H(0, 1; z)) - z'^ (8 + — H(0; + — H(l; z)) 

D o o 

17 197r2 ^ 7081 95Tr^ 32(3. ,373 Stt^ , , 38_ , , 



''7r2H(0;z) 



^H(l; 2) + iH(0, 0; ^) - ^H(0, 1; 2;) - 6H(1, 0; z) 



- 6H(1, 1; z) + 7H(0, 0, 0; z) + 6H(0, 0, 1; z) + 4H(0, 1, 0; z) + 4H(0, 1, 1; z) 
+ ]4i (f H(0; z) + ^H(0, 0; z) + ^H(0, l;z) + ^H(l, 0; 2)) + i (23 - ^ 



+ ^H(0; z) + ^H(l; z) + y H(0, 0; z) 



yH(0,l;z) + yH(l,0;z)+^H(l,l;z)) 



+ ^(y - ^ - 2C3 + ^H(0; .) - l^r^RiO; z) + ^H(l; .) + ^H(0, 0; .) 
25 

+ — H(0, 1; z) + 6H(1, 0; z) + 6H(1, 1; z) + 7H(0, 0, 0; z) + 6H(0, 0, 1; z) + 4H(0, 1, 0; z) 



+ 4H(0,l,l;z)) +z2(-— + 



32, 



g 16H(0;2;)-16H(l;z)-yH(0,0;z) 



32 32 32 

yH(0,l;^)-yH(l,0;z)-yH(l,l;z)) 



+ 0{e) 
and 



(6.10) 



+- 



+- 



e L 3 



1 + 2H(0; z) + ± + 4-1 + 2H(0; z)) - ^z^ 

^ + 5H(0; z) + 2H(1; ^) + ^ (f + ?H(0; z) + ^H(l; z)) + 6H(0, 0; 
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+ 4H(0, 1; z) + z{-^ - ^ - H(0; z) - 2H(1; z) + 6H(0, 0; z) + 4H(0, 1; z)) 



2/ 26 8, 



+ ^'(-9--3H(0;z) 



^H(l;z)) 



+ 



200 llTT^ 

L~9 6~ 



3 

- 4C3 + ^H(0; z) - iTT^RiO; z) + ^H(l; z) + 13H(0, 0; z) + 10H(0, 1; z) 



68, 

y 



20, 

y 



H(0;z) - |7r2H(0;2) - ^H(l;z) 



+ 4H(1, 0; z) + 14H(0, 0, 0; z) + 12H(0, 0, 1; z) + 8H(0, 1, 0; z) + 8H(0, 1, 1; z) 

, 1,457 IOtt^ 88 , , 88 , , 16 , ^ 16 , 
+ 4H(1, 1; + - (— - -g- + — H(0; z) + -H(l; z) + y H(0, 0; z) + y H(0, 1; z) 

32 

2 3 "^-'"^ S" -^-'-^ ^ 

- H(0, 0; z) - 2H(0, 1; z) - 4H(1, 0; z) - 4H(1, 1; z) + 14H(0, 0, 0; z) + 12H(0, 0, 1; z) 

+ 8H(0, 1, 0; z) + 8H(0, 1, 1; z) - 4(3) + ("^ + ^ - f H(0; z) - ^ H(l; z) 

- y H(0, 0; z) - y H(0, 1; z) - y H(l, 0; z) - y H(l, 1; z)) 



+ y H(l, 0; z) + y H(l, l;z))+z{-— + --- 



+ Oie). 



(6.11) 



The integrals of the one-loop virtual corrections to the {g,gq) antenna functions at 
leading colour and for a closed quark loop yield: 



V. 



9,91 



1 



+ - 



-4^(1--) 



1 r 13 c-/ 77r2 31 , , , 5,,, , 2H(0;z) 1,31 

+ 2H(0; z) + H(l; z)) + z{^- + 4H(0; z) + 2H(1; z)) + ^^(-^ - 4H(0; z) - 2H(1; z)) 



1 

+- 
e 



1 357r2 ^, 197 397r2 llCs^ ,37 77r2. , , 10^ , , ^ , , 

- ^H(0; - ^H(l; z) - 15H(0, 0; z) - y H(0, 1; - y H(l, 0; z) - 5H(1, 1; z) 

1 29 1 37 liT^ 29 

+ ^ (y H(0; z) + 6H(0, 0; + 3H(0, 1; + 3H(1, 0; z)) +-(-- — + -H(0; z) 



z^ 8 



12 



6 



+ y H(l; z) + 6H(0, 0; z) + 3H(0, 1; z) + 3H(1, 0; z) + 2H(1, 1; ^)) + ^(^ - ^ 
+ ^H(0; z) + ^H(l; z) + 12H(0, 0; z) + 6H(0, 1; z) + 6H(1, 0; z) + 4H(1, 1; z)) 



^ - y H(0; 2) - y H(l; z) - 12H(0, 0; z) - 6H(0, 1; z) - 6H(1, 0; z) - 4H(1, 1; z)) 



+ 



24l7r2 125C3 ^ 239 2837r2 IItt^ 793(3. . 



48 



22C; 



, , / 63 77r^s„ , , ^>^„ 
)^o(^) + (-y + -^)l?i(^) + Y^^: 



12 96 

29, 



144 

2 



317 _ llTT^ 



•2[z) - -V,{z) - ^H(0; z) + f vr^HCO; z) - Hh(1; ^) 
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10 19 37 37 43 

+ y 7r2H(l; z) - -H(0, 0; z) - -H(0, 1; z) - -H(l, 0; z) - -H(l, 1; z) 

25 

- 35H(0, 0, 0; z) - 20H(0, 0, 1; z) - 15H(0, 1, 0; z) - yH(0, 1, 1; z) - 15H(1, 0, 0; z) 



25, 



1 /TT 



91, 



(ir-C3 + ^H(0;z) 



- y H(l, 0, 1; z) - 10H(1, 1, 0; z) - 10H(1, 1, 1; z) + ^ _ ^ , ^ 

7 1 37 19 47 

- -7r2H(0; z) - -^'H(1; z) + -H(0, 0; z) + — H(0, 1; z) + -H(l, 0; z) + 14H(0, 0, 0; z) 
3 fa 3 3 

+ 8H(0, 0, 1; + 6H(0, 1, 0; z) + 5H(0, 1, 1; z) + 6H(1, 0, 0; z) + 5H(1, 0, 1; z) 

Ca 



+ 4H(1,1,0;2;)) +-( 



28 1Q 29 29 

_^H(0, 0; z) + ^H(0, 1; z) + -f H(l, 0; z) + -f H(l, 1; z) 

3 3 D D 



+ 14H(0, 0, 0; z) + 8H(0, 0, 1; z) + 6H(0, 1, 0; z) + 5H(0, 1, 1; z) + 6H(1, 0, 0; z) 
+ 5H(1, 0, 1; z) + 4H(1, 1, 0; z) + 4H(1, 1, 1; z)) + ^(^^ - - ^ + ^H(0; z) 
- y vr^HlO; + ^H(l; z) - z) + y H(0, 0; z) + y H(0, 1; z) + y H(l, 0; z) 

+ ^H(l, 1; z) + 28H(0, 0, 0; z) + 16H(0, 0, 1; z) + 12H(0, 1, 0; z) + 10H(0, 1, 1; z) 



a77r^ 



+ 



50C3 



+ 12H(l,0,0;z) + 10H(l,0,l;z)+8H(l,l,0;z)+8H(l,l,l;z)) +,. ^ ^ , ^ 

+ y 7r2H(0; z) + ^7r2H(l; z) - y H(0, 0; z) - y H(0, 1; z) - y H(l, 0; z) - y H(l, 1; z) 

- 28H(0, 0, 0; z) - 16H(0, 0, 1; z) - 12H(0, 1, 0; z) - 10H(0, 1, 1; z) - 12H(1, 0, 0; z) 

- 10H(1, 0, 1; z) - 8H(1, 1, 0; z) - 8H(1, 1, 1; z)) 
+ 0{e) (6.12) 



and 



+ - 



= 

9,91 




1 


- z) 




1 




-z)- ^Vo{z 


1 - 

€ - 







3i-r + r 



5.w^ X H(0;z) I, 1 1,,^ , 
6^(«'^)-3(H)+i(-4-3^(«'^) 



- iH(l; z)) + ^(-3^ - ^H(0; z) - ^H(l; z)) + z'{hi{Q- z) + ^H(l; ^)) + ^H(l; z) 



+ 



1 57r2 ^,7 37r2 

2- — + ^^^-^^6-1^ 



9 



- iH(l; z) + ^H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z) + ^H(l, 1; z) + ^ (-^H(0; z) 

- Jh(0,0;z) - Jh(0,1;z) - Jh(1,0;z)) + -(-7 + 



y-iH(0;.)-lH(l;.) 



-51 - 



- ^H(0, 0; z) - ^H(0, 1; z) - ^H(l, 0; z) - ^H(l, l-z))+z{-^ + '^- ^H(0; z) 

- ^H(l; z) - ^H(0, 0; z) - ^H(0, 1; z) - hl{l, 0; z) - ^H(l, 1; z)) + (- y 
+ ^H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z) + ^H(l, 1; z)) 



+ 0{e). 



(6.13) 



6.3 Gluon-gluon antennae 

The NLO gluon-gluon antenna functions crossed to initial-final kinematics integrate to: 



9,99 



1 

+ - 

1 

+- 
e 



2(5(1 - z) 



4+^d{l-z) 



- 2z + 2z'^ - 2Vo{z) 



+ [^(1 - ^) (l^ - y) - y ^0(2) + 2Pi(z) + 4H(0; z) + 4H(1; z) - ^H(0; z) 



+e 



- ^ (^ + 2H(0; z) + 2H(1; z)) + (-z + z^) (2H(0; z) + 2H(1; z)) 

11. 



2 - TT^ + <5(1 - .) (^ - - ^) + (-^ + ^)Po(.) + ^V,{z) - V,{z) + 4H(0, 0; .) 



18 
1 



2 

11, 



+ 4H(0, 1; z) + 4H(1, 0; z) + 4H(1, 1; 2) + - (-- g-H(0; z) - 2H(0, 0; z) - 2H(0, 1; z 
- 2H(1, 0; z)) + i (-^ + ^ - ^H(0; z) - ^H(l; z) - 2H(0, 0; z) - 2H(0, 1; z) 



- 2H(1, 0; z) - 2H(1, 1; z)) + {z - z^) (y - 2H(0, 0; z) - 2H(0, 1; z) - 2H(1, 0; z) 
-2H(l,l;z)) 



6 



28C3 



3 

,67 7r2 



+ <5(l-z)( 



1214 677r2 tt'' 77C3 



202 llvr^ 14C: 



81 



72 



48 



+ 



+ - y)^^i(-^) - ^^^2(^) + JP3(^) + 2H(0; z) - 7r2H(0; z) + 2H(1; z) - 7r2H(l; z) 

4H(0, 0, 0: z) + 4H(0, 0, 1; z) + 4H(0, 1, 0; z) + 4H(0, 1, 1; z) + 4H(1, 0, 0; z) 



1 



67, 



1 



+ 4H(1, 0, 1; z) + 4H(1, 1, 0; z) + 4H(1, 1, 1; z) + ^ (-^^^(O; z) + -7r2H(0; z 

- y H(0, 0; z) - y H(0, 1; z) - y H(l, 0; z) - 2H(0, 0, 0; z) - 2H(0, 0, 1; z) 

- 2H(0, 1, 0; z) - 2H(0, 1, 1; z) - 2H(1, 0, 0; z) - 2H(1, 0, 1; z) - 2H(1, 1, 0; z)) 

- y H(0, 0; z) - y H(0, 1; z) - ^ H(l, 0; z) - y H(l, 1; z) - 2H(0, 0, 0; z) 

- 2H(0, 0, 1; z) - 2H(0, 1, 0; z) - 2H(0, 1, 1; z) - 2H(1, 0, 0; z) - 2H(1, 0, 1; z) 



-52- 



- 2H(1, 1, 0; z) - 2H(1, 1, 1; z)) + [z - z") + ^tt^HIO; z) + ^7r2H(l; z) 

- 2H(0, 0, 0; z) - 2H(0, 0, 1; z) - 2H(0, 1, 0; z) - 2H(0, 1, 1; z) - 2H(1, 0, 0; z) 

- 2H(1, 0, 1; z) - 2H(1, 1, 0; z) - 2H(1, 1, 1; z)) 

for the (5, 55) case and to 



(6.14) 



- = 



1 

+- 
e 



+e 



<5(1 - z) 



-^.(i-.)+ipo(.)+^: 



H(0;z) + i(^ + iH(0;z) + ^H(l;^)) 



164 57r2 7C3x ,28 7r2, , , 5 , , 1 , , 1 ,5 , , 

+ ^H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z)) + ^ (P - ^ + ^H(0; z) + ^H(l; z) 
+ ^H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z) + ^H(l, 1; z)) 



(6.15) 



for the (5, qq) case. 

The integrated form of the NNLO double real radiation antenna function for {g, ggg) 

is: 



9,999 



1 

1 r 



105(1 



20 



1 



40 + 225(1 - z) - 20X'o(^) - — - 20z + 20z^ 

68 + 5(1 -z)(— -Stt^) -44X'o(^) +40r»i(z) + 64H(0;z) + 80H(l;z) - ^^R{0;z) 
9 ^ 

1 220 

- - (— + 40H(0; z) + 40H(1; z)) + z{-46 - 56H(0; z) - 40H(1; z)) 

Z o 

154 

+ (— + 40H(0; z) + 40H(1; z)) 



1 

+- 
e 



4868 1547r2 272(3 



1180 



) + ( — + l6n^)Vo{z) + 88Vi{z) - AOViiz) 



^-28^' + 6{l-z){ 
482 

+ — H(0; z) + 136H(1; z) - 16H(-1, 0; z) + 128H(0, 0; z) + 128H(0, 1; z) 

+ 160H(1, 0; z) + 160H(1, 1; z) + (-88H(0; z) - 84H(0, 0; z) - 80H(0, 1; z 
1 „ , I, 1556 467r2 



80H(1, 0; z)) + —— (— + 8H(-1, 0; z) - 4H(0, 0; z)) + -{ — + 

i. -\- Z o Z V 
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352 440 

- — H(0; z) - — H(l; z) - 8H(-1, 0; z) - 80H(0, 0; z) - 80H(0, 1; z) - 80H(1, 0; z) 

- 80H(1, 1; z)) + + ^ - ^H(0; z) - 92H(1; z) - 8H(-1, 0; z) 

- 128H(0, 0; z) - 112H(0, 1; z) - 80H(1, 0; z) - 80H(1, 1; z)) + - 

9 o 

308 

+ 132H(0; z) + — H(l; z) - 8H(-1, 0; z) + 88H(0, 0; z) + 80H(0, 1; z) + 80H(1, 0; z) 
+ 80H(1,1;2)) 

r5249 5927r2 1424(3 ,,25811 14087r2 3197r^ 934C3x _ , w 9736 
[-9 9 3- + - ")(^4 2^ + ^80 ^) + 

3087r2 544(3 X .2360 / n / n 80_ , , 32 2... , ^ 1759 , 

+ — ^ + -^) + (-^ - 327r2)l?i(^) - 88P2(-2) + y253(^) + y7r2H(-l; 2) + -^H(0; z) 

noo 200 1 030 

- 647r2H(0; z) + — H(l; z) - — vr^HCl; z) + 16H(-1, 0; z) + — ^H(0, 0; z) 

3 3 3 

+ ^H(0, 1; z) + ^H(l, 0; z) + 272H(1, 1; z) + 64H(-1, -1, 0; z) - 80H(-1, 0, 0; z) 

- 32H(-1, 0, 1; z) - 96H(0, -1, 0; z) + 256H(0, 0, 0; z) + 256H(0, 0, 1; z) 

+ 224H(0, 1, 0; z) + 256H(0, 1, 1; z) + 272H(1, 0, 0; z) + 320H(1, 0, 1; z) + 304H(1, 1, 0; z) 

+ 320H(1, 1, 1; + ^ (-^ + 44C3 - ^H(0; z) + 407r2H(0; + ^vr2H(l; z) 

- — H(0, 0; z) - 176H(0, 1; z) - — H(l, 0; z) + 24H(0, -1, 0; z) - 172H(0, 0, 0; z) 

- 168H(0, 0, 1; z) - 144H(0, 1, 0; z) - 160H(0, 1, 1; z) - 136H(1, 0, 0; z) - 160H(1, 0, 1; z) 

- 152H(1, 1, 0; z)) + ^ (+28C3 - y vr2H(-l; z) + ^tt^HCO; z) - 32H(-1, -1, 0; z) 

+ 40H(-1, 0, 0; z) + 16H(-1, 0, 1; z) + 24H(0, -1,0; z) - 12H(0, 0, 0; z) - 8H(0, 0, 1; z)) 

I, 2102 5507r2 592C3 16 o.w , x 1216„,^ , 116 n.w^ x 
+ -{ 3- + -9- + + y ^'H(-l; z) - ^H(0; z) + —irM^; z) 

- ^H(l; .) + l^vr^Hll; .) + f H(-l, 0; .) - I^H(0, 0; .) - I^H(0, 1; .) 

- — H(l, 0; z) - — H(l, 1; z) + 32H(-1, -1, 0; z) - 40H(-1, 0, 0; z) - 16H(-1, 0, 1; z) 

- 160H(0, 0, 0; z) - 160H(0, 0, 1; z) - 144H(0, 1, 0; z) - 160H(0, 1, 1; z) - 136H(1, 0, 0; z) 

- 160H(1, 0, 1; z) - 152H(1, 1, 0; z) - 160H(1, 1, 1; .)) + ^ 31br^ ^ 688^ 
+ f -^H(-l; .) - ^H(0; z) + ip;r^H(0; .) - ^H(l; z) + H,2h(i. ,) 

+ 16H(-1, 0; - — H(0, 0; z) — H(0, 1; z) - — H(l, 0; z) - 184H(1, 1; z) 

+ 32H(-1, -1, 0; z) - 40H(-1, 0, 0; z) - 16H(-1, 0, 1; z) - 272H(0, 0, 0; z) 

- 256H(0, 0, 1; z) - 208H(0, 1, 0; z) - 224H(0, 1, 1; z) - 136H(1, 0, 0; z) - 160H(1, 0, 1; z) 

- 152H(1. 1,0; .) - 160H(1, M; .)) + .^(^ - - ^ + f .=H(-l; .) 
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+ M?H(0; .) - f .=H(0; .) + if H(l; .) - M..„(i; .) + |h(-1, 0; .) 

880 , , , , 616 , , 616 , 

+ — H(0, 0; z) + 264H(0, 1; z) + — H(l, 0; z) + — H(l, 1; z) 

+ 32H(-1, -1, 0; z) - 40H(-1, 0, 0; z) - 16H(-1, 0, 1; z) - 48H(0, -1, 0; z) 

+ 184H(0, 0, 0; z) + 176H(0, 0, 1; z) + 144H(0, 1, 0; z) + 160H(0, 1, 1; z) + 136H(1, 0, 0; z) 

+ 160H(1, 0, 1; z) + 152H(1, 1, 0; z) + 160H(1, 1, 1; z)) 

+ 0{e). (6.16) 



The [g, qqg) antenna function has a leading and subleading colour contribution. They 
integrate to: 



-.0 

^gmg 



+- 



+ - 



+- 



1 



-3^(1-^) 

7 47,,, , 10 



Vo{z) + H(0; z) + -^+z{^ + H(0; z)) - \z} 



+ 



59 ^. 55 7 237r\ 94_ , , 20_ , , 35„,^ , , 

+ + + ^) + ^^(^5 + ^) + ^) + ~ y 

+ ^H(0; + 3H(1; z) + 3H(0, 0; z) + 2H(0, 1; z)) + (-^ - ^H(0; z) - ^ H(l; z)) 
259 llTT^ ^1 28613 257r2 IISCsn /55 7 237r2. , , 188^ 



12 4 ^ 648 ■ 6 

20_ , , 113„,^ , 7 , 59 

+ -P2(^) - -^H(0; z) - -7r2H(0; z) ^ - 



9 ' ■ M8 9 

59 35 
H(l;z)--H(0,0;z)--H(0,l;2;) 
D 3 



40 1 ItP' 188 44 40 

- y H(l, 0; z) - 14H(1, 1; + ^ (— + — H(0; z) + y H(0, 0; z) + y H(0, 1; z) 

44,w, „ 1.3065 lOvr^ 296,,,^ , 296^^,, , ..tt/^ n x 
+ y H(l, 0; z)) + -{—^ — + -^H(0; z) + ^H(l; z) + 16H(0, 0; z) 

+ 16H(0, 1; z) + 16H(1, 0; z) + 16H(1, 1; z)) + 7H(0, 0, 0; z) + 6H(0, 0, 1; z) 

+ 4H(0, 1, 0; z) + 4H(0, 1, 1; ^) + - ^ - 2(3 + ^H(0; z) - ^7r2H(0; z) 



41 49 
+ -H(l;z) + -H(0,0;z) 



y H(0, 1; z) + y H(l, 0; z) + 6H(1, 1; z) + 7H(0, 0, 0; z) 

56, 



757 207r^ 

+ 6H(0, 0, 1; z) + 4H(0, 1, 0; z) + 4H(0, 1, 1; z)) + z^ (__ + __ _ ^ H(0; z) 

56 32 39 32 32 

- y H(l; z) - y H(0, 0; z) - y H(0, 1; z) - y H(l, 0; z) - y H(l, 1; z)) 



(6.17) 



and 



1 








1 




+- 


1 


1 r 


7 


+- 




e L 


3 " 



19 4 8 4 

1 - -5(1 -z) + -Voiz) + 2H(0; z) + - + z{-l + 2H(0; z)) - - 



3 



373 57r^ 



38. 
9 ' 



(z) - -Vi{z) + 3H(0; z) + 2H(1; z) + 6H(0, 0; z) 



+ ^) + Irrf + i + y ^^(0^ ^) + y H(i; .)) + _ ^ - h(o; .) 

- 2H(1; z) + 6H(0, 0; z) + 4H(0, 1; z)) + ^^(-^ _ ^H(0; z) - ^H(l; z)) 
391 137r2 ^, 6973 95^2 54^ 373 10^2 75 



8 52 7 14 25 

+ -I?2(^) + y H(0; z) - -7r''B{0; z) + y H(l; z) + y H(0, 0; z) + 6H(0, 1; ^) 

16 

+ — H(l, 0; z) + 4H(1, 1; z) + 14H(0, 0, 0; z) + 12H(0, 0, 1; z) + 8H(0, 1, 0; z) 

O 

1 7G 16 16 

+ 8H(0, 1, 1; z) + ^ (-H(0; z) + 8H(0, 0; z) + y H(0, 1; z) + y H(l, 0; z)) 

1 ,809 287r2 152^^, , 152^^, , 32^^, 
+ - - + -^H(0; z) + — H(l; z) + y H(0, 0; z) 

+ y H(0, 1; z) + 8H(1, 0; z) + y H(l, 1; z)) + z{-^ _ l!^ _ 4(3 + ^H(0; z) 



- j7r2H(0; - ^H(l; 2) - ^H(0, 0; z) - 2H(0, 1; z) - ^H(l, 0; z) - 4H(1, 1; z) 



14, 



11, 



20, 



+ 14H(0, 0, 0; z) + 12H(0, 0, 1; z) + 8H(0, 1, 0; z) + 8H(0, 1, 1; z)) + (-^ + ^ 



-fH(0;.) 
+ 0(6). 



y H(l; z) - y H(0, 0; - y H(0, 1; - y H(l, 0; z) 



16 



H(l,l;z)) 



(6.18) 



Finally, the integrated one-loop antenna functions for (5, gg) at leading colour and for 
a closed quark loop read: 



9,99 



+ 



1 
1 



5{l - z) 



55 1 
-2 - —5(1 -z)+ Vq{z) + - + z-z'^ 



+ ^ (y + 4H(0; + 2H(1; z)) + -^(y + 4H(0; z) + 2H(1; z)) 



+ z' 



1 

+- 
e 



11 

~Y 

77r2 



4H(0;z) -2H(l;z)) 



c-/, X/ 196 7 3857r2 22(3^ ,389 77r2, ^ ^ 22_ , , , , 

22 22 

- — H(0; z) - — H(l; - 24H(0, 0; z) - 12H(0, 1; - 12H(1, 0; z) - 8H(1, 1; z) 
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+ 



11, 



(11H(0; z) + 12H(0, 0; z) + 6H(0, 1; z) + 6H(1, Q- z)) + z[^ - ^ + ^H(0; z) 

L — z boo 



11, 



77r^ 



+ 



+ — H(l; z) + 12H(0, 0; z) + 6H(0, 1; z) + 6H(1, 0; z) + 4H(1, 1; 2)) +z'^{ „ 

O 

- y H(0; z) - y H(l; z) - 12H(0, 0; z) - 6H(0, 1; z) - 6H(1, 0; z) - 4H(1, 1; z)) 

1 389 Ttt^ 22 

+ - (-3g- - -g- + 11H(0; z) + y H(l; + 12H(0, 0; z) + 6H(0, 1; ^) + 6H(1, 0; z) 

+ 4H(l,l;z)) 

1697r2 lOOC: 



-13 



18 



+ 



3 14453 31557r2 IItt^ 269(3^ ,3197 557r2 



324 ' 432 
11 



^)Vo{z) + (-^ + ^)Pi(z) + ^P2(^) - ^Vs{z) - ^H(0;z) + ^7r2H(0; z) - 4H(l;z) 



108 
28 

T 



+ y vr2H(l; z) - y H(0, 0; z) - |h(0, 1; z) - y H(l, 0; z) - |h(1, 1; z) 

- 56H(0, 0, 0; z) - 32H(0, 0, 1; 2) - 24H(0, 1, 0; z) - 20H(0, 1, 1; z) - 24H(1, 0, 0; z) 

- 20H(1, 0, 1; z) - 16H(1, 1, 0; z) - 16H(1, 1, 1; z) + - 2(3 + ^H(0; z) 

14 1 88 44 5^ 

- -TT^m z) - -7r2H(l; z) + y H(0, 0; z) + y H(0, 1; z) + y H(l, 0; z) 

+ 28H(0, 0, 0; z) + 16H(0, 0, 1; 2) + 12H(0, 1, 0; z) + 10H(0, 1, 1; z) + 12H(1, 0, 0; z) 
+ 10H(1. 0, l; .) + 8H(1. 1, 0; .)) + i - - + ?|h(0; .) 

- ^7r2H(0; z) + ^H(l; z) - ^vr2H(l; z) + 22H(0, 0; z) + ^H(0, 1; z) + 11H(1, 0; z) 
+ 11H(1, 1; z) + 28H(0, 0, 0; z) + 16H(0, 0, 1; z) + 12H(0, 1, 0; z) + 10H(0, 1, 1; z) 

+ 12H(1, 0, 0; z) + 10H(1, 0, 1; + 8H(1, 1, 0; z) + 8H(1, 1, 1; ^)) 

- ^ + ^H(0; - y TT^HIO; z) + ^H(l; z) - ^tt^HII; z) + ^H(0, 0; z) + y H(0, 1; z) 

+ Jh(1, 0; z) + ^H(l, 1; z) + 28H(0, 0, 0; z) + 16H(0, 0, 1; z) + 12H(0, 1, 0; z) 
+ 10H(0, 1, 1; z) + 12H(1, 0, 0; z) + 10H(1, 0, 1; z) + 8H(1, 1, 0; z) + 8H(1, 1, 1; z)) 

+ + ^ + y vr^HCO; .) + ^vr^Hll; .) - Hh(0, 0; .) - Hh(0, 1; .) 

- ^H(l, 0; z) - ^H(l, 1; z) - 28H(0, 0, 0; z) - 16H(0, 0, 1; z) - 12H(0, 1, 0; z) 



- 10H(0, 1, 1; z) - 12H(1, 0, 0; z) - 10H(1, 0, 1; z) - 8H(1, 1, 0; z) - 8H(1, 1, 1; z)) 
+ 0{e) 
and 
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+ - 
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6{l-z){ 



143 57r\ 11 



2H(0;2) 
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2^8 



+ ^H(0; z) + ^H(l; z)) - z(i + ^H(0; z) + ^H(l; z)) + (^H(0; z) + ^H(l; z)) 



4 57r2 ^,979 SStt^ 



14C3. . 38 57r2 , . 11^ . . 1^ . . 



+ ^H(0, 0; z) + ^H(0, l-z) + ^H(l, 0; z) + ^H(l, 1; z) + ^ (" ^^^(0; z) 
- ^H(0, 0; z) - |h(0, 1; - f H(l, 0; z)) + i (-g + ^ - iiH(0; z) - ^H(l; z) 
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- ^H(0, 0; z) - ^H(0, 1; - |h(1, 0; z) - ^H(l, 1; z)) + , 



18 

7 57r2 

'9 6 



1 



H(0;z) 



1 2 2 2 2 57r^ 

- -H(l;^) - -H(0,0;^) - -H(0,1;^) - -H(l,0;z) - -H(l, 1; ^)) + (- — 

+ ^H(0, 0; z) + ^H(0, 1; z) + ^H(l, 0; z) + ^H(l, 1; z)) 

+ Oie) , 



(6.20) 



while the leading and subleading colour and quark loop contributions to the one-loop (5, qq) 
antenna function integrate to: 
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+- 



+ 
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5(1 - z) 



-18 

119,,, X 11^ / X 

<^(1 - ^) - TZ^oiz) 
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2H(0;z) I, 11 2^^, 



9 +^<l-'-H324 



787 557r2. 



32, 



11 



-^2^o(^) + Y^^?i(^)- 3H(0;^) 



1 
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- ^H(l,0;z) + ^(-!^ - §H(0;z) - ^H(0,0;^) - ^H(0,l;z) - ^H(l,0;z)) 



H(0,0;z) 



1-z 

+ 7 (-S + ^ - ^) - T^H(1; - 2H(0, 0; z) - ^H(0, 1; z) - ^H(l, 0; z)) 



TT 

108 ' "9 

,27r2 2, 



18 



+ + -H(0,0;2:) + -H(l,0;z)) 
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and 
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1 

+- 
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5^(1-.) 

19 

— 5(l-z) 
36 ^ ' 



3 ^ 3z. 



,,173 57r2, 19^ , , 2^ , 

5 1-2)( ) Vq{z) + -Vi{z, 

^ ^^408 36 ^ 18 ^ ^ 3 ^ ' 3(1 



2H(0:z) I, 19 2^^,^ , 2„,, 



-58- 



+ 



2 ^2 343 957r2 TCsx / 173 Svr^ ^ ^ 19^ ^ ^ 2^ , , 



-^H(0; 0)-^H(O,O;3) 



1 1 19 4 

-H(l, 0; z) + ^ (-- g-H(0; z) - 2H(0, 0; z) - -H(0, 1; z) 



- ^H(l, 0; z)) + i (-^ + ^ - ^H(0; z) - ^H(l; z) - ^H(0, 0; z) - ^H(0, 1; z) 



2 4 , ,27r2 2 2 

-H(l, 0; - -H(l, 1; z)) + z{— + -H(0, 0; z) + -H(l, 0; z)) 



and 
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5{l-z)--V^{z) 
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,,53 Stt^ 5_ , , 1_ . . H(0;^) 1, 5 1„,^ , 1^^, 
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7. Rederivation of NNLO coefficient functions 

Being derived from physical matrix elements, the integrated antenna functions can be 
compared to results from literature for inclusive cross sections or coefficient functions, 
as was done previously for the final-final antennae in [25,38]. In the case of the initial- 
final antennae, we can compare the integrated quark-antiquark antennae and gluon-gluon 
antennae against NNLO corrections to deep inelastic coefficient functions known in the 
literature. The former ones can be checked against DIS structure function calculations [51] 
whereas the latter can be compared to the <^DIS structure functions computed in [55,56]. 
The quark-gluon antennae, derived from neutralino decay, cannot be associated to any 
physical process and only the deepest pole could be checked against a combination of 
Altarelli-Parisi splitting functions. 

7.1 Deep inelastic scattering: 7 induced 

In deep inelastic scattering (DIS) the hadronic tensor W^'^ is contracted with the metric 
tensor — 5^,^. This corresponds to the trace of the hadronic tensor, which in terms of the 
structure functions J^l and is given by 

-Wi: = -^^L {z, Q') + ^:F, {z, Q') , (7.1) 

where the structure functions can be expanded in powers of the strong coupling constant, 
following the notation of [51] (but using d = 4 — 2e, as throughout the rest of this work): 
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and = cxp [e (In (47r) —7b)]. For clarity, we drop the dependence on z and Q^. The 
factors Se and fi are the conventional factors appearing in dimensional regularization. 

To zeroth order in the structure functions are given by the simple parton model 
result 

^^lI = ^^lI = ^^ J-(2=^(1-z), -^S=0. (7.3) 

Since the overall normalization of the antenna functions is given by the leading order 
antenna Aq^q = (5(1 — z), we find that the correct normalization of Wji to be checked 
against the antennae is 

2 WI: = J^2-^J'l. (7.4) 



d-2 " d-2' 

This last equation at order ag and can be compared to a linear combinations of NLO 
and NNLO antennae respectively. For completeness and clarity, before giving the explicit 
relations between the structure functions and the antennae, we report the expressions of 
the structure functions in terms of the Altarelli-Parisi splitting functions and the other 
coefficient functions Cj and Oj. We take them from [51], adjusting to d = 4 — 2e and 



(1) 



^..-^K-S-^2a«)' (7-5) 

4'i=^i^(-^n^?+4:]-^2ag), (7.6) 

4^(ca-2aa), (7.7) 

4'i = (-^<^ + 4;^^2a«), (7.8) 

+ - 4/?oaa + P£) a« + P^f a^) , (7.9) 

+ eg - 4/3oa« + Pi°) + P.^ 4;]) , (7.10) 

= (7 (^/^ocS - ^ oil) + cgf - 4/3oaa + P(0) ^ a« ) , 

(7.11) 



(2)NS _ (^1 f 1 p(0) ^ p(0) _ ^ piO)\ 
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(2)PS 



L,q 
(2)PS 



2,? 



-2,? 



2e 
1 /I 



9ff 



,(1) 
2,ff 



(7.12) 
(7.13) 

(7.14) 



The previous expressions are needed to check ah the initial-final quark-antiquark antennae. 
They are still unrenormalized and are thus to be compared to the unrenormalized antenna 
functions. The following identities hold at NLO: 



r(l) 



3-26_^(i) 



Cf V ^'^ 2 - 2e 
(2-26) f^a) _ 3-26^^(1) 



:4^,,,, + 85(1-z)fW|^ , 



Nf 



2,ff 2 - 2e ^'^ 



9,QQ ' 



(7.15) 
(7.16) 



where Fq^^ is the one-loop quark form- factor given in [57]. The notation means that 
only the terms of X proportional to a are considered. The corresponding identities at 
NNLO are 
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<2,NS) 



3-2e 
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(2,NS) 
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2,9 
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JV 

2Alggq + 8^^;,^, + 165 (1 - z) (pf + CfF^'^ 
2^0 ,,,,, + 16,5 (1 - z) (Pf + 
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(2 - 2e) 
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(2-26) 
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-2e 
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- 2e 



-(2,PS) 



,(2) 



(2) 



AT 



- ^•"9,9'9'9 ' 



4^0 _ - 8^^'^- 



(7.19) 
(7.20) 

(7.21) 
(7.22) 



The structure functions (7.5)-(7.14) are given in [51] in unrenormalized form. To com- 
pare with our results, we thus considered the one- loop unrenormalized antenna functions 
(obtained by undoing the renormalization) the two-loop unrenormalized form factor, both 
characterized by the label U. Full agreement is found with [51]. 
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7.2 Deep inelastic scattering: 0-scalar induced 

The gluon-gluon T-^ Q- and ?Y-type antenna functions can be checked with the deep in- 
clastic scalar induced structure functions derived in [55,56]. These structure functions are 
obtained in an effective theory with a scalar ^ coupled to the square of the gluon field 
strength tensor. This effective theory is very instrumental in deriving gluonic splitting 
functions, and is also the effective field theory describing the coupling of the Higgs boson 
to gluons through a loop of asymptotically heavy coloured particles. Since is a scalar 
particle the corresponding hadronic tensor has only one structure function which can be 
either quark- or gluon-initiated. Their perturbative expansion reads: 




(7.23) 



At NLO and NNLO the structure functions have the following form: 



^(2) -(^{\ p(0) ^ p(0) , ^ p(0) ^ p(0) , ^ p(0)\ 
^4>,9 IS^^ ^^99 ^ g -^99 ^ ^99 ^ 2 ^99 J 

+ 1 (ApW _ i p(0) ^ c^i) _ ^ p(o) ^ g( 

The explicit expressions for the coefficients c^^], a^^] and c^^], with i = q,g are given in 
Appendix B. They were derived independently in [56], and wc find full agreement. In 
our calculation, they follow from the following identities between structure functions and 
antennae: 

= 2 (C^^,„ + JVFS»,f) + « (1 - .) Fl') , (7.28) 
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r(2) 



N2 



(2) 



(2) 



1 .^(2) 



7V2 



^ - Qlm + + 45 (1 - .) (2Ff ) + ^) 



Ar2 



1 -j(2) 



N 



1 



(2) 



(7.30) 

(7.31) 
(7.32) 
(7.33) 

(7.34) 

(7.35) 

(7.36) 



The antenna Jq^qqq comes from the interference of four quark final states with identical 
quark flavour, which is not contained in Table 1 since it is finite in all limits. It is given by 

^Im = (-i + i + I + (2 - ^ - -)C3) + O (^) • (7.37) 

The functions F^*'* (i = 1,2) are the gluon form factor at one- respectively two- loop given 
in [57]. 



8. Conclusions and Outlook 

In this paper, we have extended the NNLO antenna subtraction formalism [24] to include 
initial-final antenna configurations, where one of the hard radiator partons is in the initial 
state. We described the construction of initial-final antenna subtraction terms from antenna 
functions, including the required phase space factorization and mappings. 

The analytic integration of the NNLO initial-final antenna functions requires the 
derivation of inclusive 2^3 phase space integrals at tree-level and of 2 ^ 2 phase space 
integrals at one loop. We reduced all these integrals to a small set of master integrals, which 
we computed using differential equations or by direct integration. With these, we provide 
integrals of all NNLO initial-final antenna functions, as required for the implementation 
into a parton-level event generator programme. A highly non-trivial check of our results 
is provided by the rederivation of deep inelastic coefficient functions [51,55,56], which can 
be expressed as particular linear combinations of antenna functions. 

The subtraction terms presented here allow the construction of a parton-level event 
generator program for the calculation of NNLO corrections to jet production observables 
in deeply inelastic electron-proton scattering. Of foremost importance in this context is 
the dijet production cross section in deep inelastic scattering [58], which has been used at 
HERA to measure the gluon distribution in the proton and the strong coupling constant 
Us- This cross section is at present only known to NLO [53], which is insufficient for various 
reasons. The error on the extraction of ag from HERA (2+1) jet data is dominated [59] 
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not by the statistical uncertainty on the data, but by the uncertainty inherent to the NLO 
calculation, as estimated by varying renormalization and factorization scales. Given that 
the statistical precision of the data will further improve once all data from the HERA-II 
run are analysed, the theoretical description requires the inclusion of the next-to-next- 
to-leading order (NNLO), i.e. 0{al), corrections. Moreover, present determinations of 
parton distribution functions [60] at NNLO accuracy [61] can only include data sets for 
observables where NNLO corrections are known [51,62]. Consequently, the precise HERA 
deep inelastic dijet data could be used in global NNLO determinations only once the NNLO 
QCD corrections to this observable are computed. Using the subtraction formalism derived 
here, and the already known matrix elements [3,54,63,64], such a calculation can now be 
envisaged. 

Moreover, the initial-final antenna functions derived here are an important ingredient 
to the calculation of NNLO corrections to jet observables at hadron colliders. In this case, 
all radiator configurations contribute: final- final (derived in [24]), initial-final (derived here) 
and initial-initial. The initial- initial antenna functions are at present known only to NLO, 
and work on their extension to NNLO is ongoing. 
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A. Master Integrals 

The common prefactor Sr is defined as: 



A.l Master integrals for double real radiation 

The double real master integrals expanded to the needed order in e are given by: 




(A.l) 




(A.2) 
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2 H(0; x) - e (^TT^ - 10 H(0; x) - 6 H(0, 0; x) + 4 H(l, 0; x)) - tt^ + 4(3 

-38H(0;x) +27r^H(0;x) - ^Tr'^}l{l; x) -30H(0,0;x) +20H(l,0;a;) +4H(0,l,0;x) 

-14H(0,0,0;x) + 12H(l,0,0;x) -8H(l,l,0;x)) + (-yTr^ - 2OC3 + ^tt^ 

20 

+ 130 H(0; x) - IOtt^ H(0; x) - 24 H(0; x) (3 + ^tt^ H(1; x) + 8 H(1; x) Cs 

+ ^7r2H(0,l;x) + 114H(0,0;x) - ^7r2H(0,0;x) - 76H(1, 0; x) + 47r2 H(l, 0; x) 

- ^7r2H(l,l;x) -20H(0,l,0;x) - 4H(0,0, l,0;x) + 70H(0,0,0;x) - 60H(l,0,0;x) 

+ 40H(l,l,0;x) +8H(l,0,l,0;x) - 12 H(0, 1, 0, 0; x) + 8H(0, 1, 1, 0; x) 

+ 30H(0,0,0,0;x) - 28H(1, 0, 0, 0; x) + 24H(1, 1, 0, 0; x) - 16H(1, 1, 1, 0; x)) 

+ O (6^)] , 

(A.3) 

l26 =Sr{QT''^'il-x)-^^x7r 

+ 2 H(0, 0; x) + 2 H(l, 0; x)) + 4(3 - ^tt^ H(0; x) - ^tt^ H(1; x) 
-4H(0,l,0;x) +6H(0,0,0;x) +6H(l,0,0;x) -4H(l,l,0;x) +e(^7r^ 
- 4H(0; x)C3 - 4H(1; x)C3 + ^tt^ H(0, 1; x) - 27r^ H(0, 0; x) - 2Tr^ H(l, 0; x) (^-4) 

+ ^TT^ H(l, 1; x) - 4H(0, 0, 1, 0; x) - 4H(1, 0, 1, 0; x) - 12H(0, 1, 0, 0; x) 

+ 8H(0, 1, 1, 0; x) + 14H(0, 0, 0, 0; x) + 14H(1, 0, 0, 0; x) - 12H(1, 1, 0, 0; x) 
+ 8H(l,l,l,0;x)) +0(e^)] , 

/249 =Sr {qy-' {1 - xr''x\ 

2 8 1 / 10 \ 

— 3 - ^ H(0; x) + - ( — TT^ + 10 H(l, 0; x) - 22 H(0, 0; x) ) + 56(3 + IOtt^ H(0; x) 

6 6 6 V O / 

- ^7r2H(l; x) + 20 H(0, 1, 0; x) - 50 H(0, 0, 0; x) + 38 H(l, 0, 0; x) - 20 H(l, 1, 0; x) 

/ 56 20 
+ 4-77^^ + 156C3H(0;x) -68H(l;x)C3 - — tt^ H(0, 1; x) + 227r2 H(0, 0; x) 
V 45 6 

20 

- 187r2 H(l, 0; x) + — tt^ H(1, 1; x) + 28 H(0, 0, 1, 0; x) - 52 H(l, 0, 1, 0; x) 

o 

+ 76H(0,l,0,0;x) - 40H(0, 1, 1, 0; x) - 106 H(0, 0, 0, 0; x) + 94 H(l, 0,0,0; x) 

- 76H(l,l,0,0;x) + 40 H(l, 1, 1, 0; x)) + O (e^) , 

(A.5) 
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/i25=-5r(Q') ' '^(l-x)-2^(l + x)-V7r 

4h(0; x) + - (^-tt^ - 4H(-1, 0; x) + 8H(0, 0; x) - 4H(1, 0; x)) - 14(3 
o 

- -TT^ H(0; x) + 27r2 H(l; x) - 8H(-2, 0; x) + 27r2 H(-l; x) - 8H(0, 1, 0; x) 
+ 8H(-1, -1, 0; x) - 16H(-1, 0, 0; x) + 8H(-1, 1, 0; x) + 20H(0, 0, 0; x) 

+ 8H(l,-l,0;x) - 16H(l,0,0;x) +8H(l,l,0;x) +e( + 47r2H(0, -1; x) 

V 30 

16 

+ 47r^H(0, 1; x) - 47r^H(-l, -1; x) + —7r'^B.{-l, 0; x) - i7r'^B{-l, 1; x) 

O 

- ^TT^RiO, 0; x) - 47r^H(l, -1; x) + ^7r2H(l, 0; x) - 47r2H(l, 1; x) 

+ 28C3H(-1; x) - 28C3H(0; x) + 28C3H(1; x) - 32H(0, 0, -1, 0; x) 

- 8H(0, 0, 1, 0; x) + 16H(0, -1,-1, 0; x) - 32H(0, -1, 0, 0; x) + 16H(0, -1, 1, 0; x) 
+ 16H(-1, 0, -1, 0; x) + 16H(-1, 0, 1, 0; x) + 16H(1, 0, -1, 0; x) + 16H(1, 0, 1, 0; x) 
+ 16H(0, 1, -1, 0; x) - 32H(0, 1, 0, 0; x) + 16H(0, 1, 1, 0; x) - 16H(-1, -1,-1, 0; x) 
+ 32H(-1, -1, 0, 0; x) - 16H(-1, -1, 1, 0; x) - 40H(-1, 0, 0, 0; x) 

- 16H(-1, 1, -1, 0; x) + 32H(-1, 1, 0, 0; x) - 16H(-1, 1, 1, 0; x) + 44H(0, 0, 0, 0; x) 

- 16H(1, -1, -1, 0; x) + 32H(1, -1, 0, 0; x) - 16H(1, -1, 1, 0; x) - 40H(1, 0, 0, 0; x) 

- 16H(1, 1, -1, 0; x) + 32H(1, 1, 0, 0; x) - 16H(1, 1, 1, 0; x)) + O (e^) , 

(A.6) 
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X TT 



H((|x) _ 7H(0,0;x) ^ ^ 7 o _ ^^^^^ _^ ^ ^^^^^^ q. 

6 6 o 

- 19H(0, 0, 0; x) + e + 20C3H(0; x) - 32H(0, 0, -1, 0; x) + 4:Tr^}i{0, -1; x) 
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- 2tt^R{0, 1; x) + ^7r2H(0, 0; x) + 10H(0, 0, 1, 0; x) + 16H(0, -1, -1, 0; x) 

- 32H(0, -1, 0, 0; x) + 16H(0, -1, 1, 0; x) + 16H(0, 1, -1, 0; x) + 22H(0, 1, 0, 0; x) 

- 20H(0, 1, 1, 0; x) - 43H(0, 0, 0, 0; x)) + O (e^) 



(A.7) 
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/l245 =Sr ^ ''(I- + xy'x^TT 



8H(0;x) ^ l^ lOvr^ ^ 

28 . 20 



eV 3 



8H(-1, 0; x) - 28H(0, 0; x) + 16H(1, 0; x)j + 28(3 



- -TT^ H(-l; x) + —TT^ H(0; x) - — vr^ H(l; x) + 8H(0, -1, 0; x) + 24H(0, 1, 0; x) 

o o o 

+ 24H(-1, 0, 0; x) - 16H(-1, 1, 0; x) - 68H(0, 0, 0; x) - 16H(1, -1, 0; x) 

497r^ 

+ 88H(0; x)C3 - 56H(1; x)C3 

- ^7r2H(0, -1; x) - ^7r2H(0, 1; x) + 16H(0, 0, -1, 0; x) - -Tr'^}l{-1, -1; x) 



+ 56H(1, 0, 0; x) - 32H(1, 1, 0; x) + e 
i7r2H(0,-l;x)-| 



1 flQ 1 a 

- 87r2H(-l, 0; x) + —Tr^U{-l, 1; x) + —Tr'^'H.{0, 0; x) + — 7r2H(l, -1; x) 

o o o 

- —TT^B{1, 0; x) + — 7r2H(l, 1; x) + 24H(0, 0, 1, 0; x) + 16H(0, -1, -1, 0; x) 
3 3 

+ 16H(0, -1, 0, 0; x) - 16H(0, -1, 1, 0; x) + 16H(-1, 0, -1, 0; x) - 16H(-1, 0, 1, 0; x) 

- 16H(1, 0, -1, 0; x) - 48H(1, 0, 1, 0; x) - 16H(0, 1, -1, 0; x) + 80H(0, 1, 0, 0; x) 

- 48H(0, 1, 1, 0; x) - 32H(-1, -1, -1,0; x) + 16H(-1, -1, 0, 0; x) 

+ 56H(-1, 0, 0, 0; x) - 48H(-1, 1, 0, 0; x) + 32H(-1, 1, 1, 0; x) - 148H(0, 0, 0, 0; x) 

- 48H(1, -1, 0, 0; x) + 32H(1, -1, 1,0; x) + 136H(1, 0, 0, 0; x) + 32H(1, 1, -1, 0; x) 

- 112H(1, 1, 0, 0; x) + 64H(1, 1, 1, 0; x)) + O (e^) , 

(A.8) 



/2356 =5r(g') ' ''{I-X) 



-i-2^xV 



L e 



2 _ 8H(0;x) 

"3 l2 



1 / \ 28 

+ - (27r2 - 28H(0, 0; x) + 4H(1, 0; x)) + 32(3 + — tt^ H(0; x) 

-8H(0,-l,0;x) +24H(0,l,0;x) -68H(0,0,0;x) +28H(l,0,0;x) +e(|256C3 

- 1^ + 47r2H(0, -1; x) - 47r2H(0, 1; x) + —7r^R{0, 0; x) - ^Tr^B.{l, 0; x) 

18 3 3 

+ 88C3H(0; x) - 32C3H(1; x) - 16H(0, 0, -1, 0; x) + 24H(0, 0, 1,0; x) 

+ 16H(0, -1, -1, 0; x) - 32H(0, -1, 0, 0; x) + 16H(0, -1, 1, 0; x) + 32H(1, 0, -1, 0; x) 

- 40H(1, 0, 1, 0; x) + 16H(0, 1, -1, 0; x) + 96H(0, 1, 0, 0; x) - 32H(0, 1, 1, 0; x) 

- 148H(0, 0, 0, 0; x) + 76H(1, 0, 0, 0; x)) + O (e^) , 

(A.9) 



/1346 =Sr (Q') ' (1 - x) 



-'-^'x\ 



_ 10 _ 20H^ + 1 _ 40H(0, 0; x)) + 92(3 + '^n^ H(0; x) - 80H(0, 0, 0; x) 

+ e(- Jtt^ + 184C3H(0; x) + ^Tr'^U{0, 0; x) - 160H(0, 0, 0, 0; x)) + O (e^) 



(A.IO) 
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A. 2 Master integrals for one-loop single real radiation 

The virtual master integrals are given by: 



Vii =2m St {Q'^Y {I - x)-^ 

^ + 4 + H(0; + e (12 - y + 4 H(0; x) + H(0, 0; x)) + (32 - ^ - 4(3 
+ 12 H(0; x) - ^vr^ H(0; x) + 4 H(0, 0; x) + H(0, 0, 0; x)) + e'^ (^80 - 27r^ - I6C3 

- ^ + 32 H(0; x) - ^TT^ H(0; x) - 4 H(0; x) Ca + 12 H(0, 0; x) - ^tt^ H(0, 0; x) 
+ 4H(0,0,0;x) +H(0, 0,0,0; x)) +C»(e^)] , 

(A.ll) 



V2A =Vu{l - x)-' 

l + eH(0;x) + 6^ ( H(0,0;x) 



vr 
~2 



+ € 



2 

'(h(0,0,0;x) -yH(0;x)) 



(A.12) 



- + 5 + 2 H(0; x) + e (^19 - y + 10 H(0; x) + 4 H(0, 0; x) J + (^65 - — - 8C3 



+ 38 H(0; x) - -TT^ H(0; x) + 20 H(0, 0; x) + 8 H(0, 0, 0; x) + e'^ 211 - — tt^ - 4OC3 



19 



— TT^ + 130 H(0; x) - — TT^ H(0; x) - 16 H(0; x) C3 + 76 H(0, 0; x) - -tt^ H(0, 0; x) 
45 3 3 



+ 40 H(0, 0, 0; x) + 16 H(0, 0, 0, 0; x)) + O (e^) 



(A.13) 



F134 =27ri5r (Q^)"^"^' (1 - x)-^x 

^ (y + H(0, 0; x) + H(l, 0; x)) - C3 - ^tt^ H(0; x) + ^tt^ H(1; x) + 3 H(0, 0, 0; x) 
+ 4 H(l, 0, 0; x) + H(l, 1, 0; x) + e(^7r^ - 4 H(0; x) C3 - 5 H(l; x) C3 
- ^TT^ H(0, 0; x) - ^TT^ H(l, 0; x) + ^tt^ H(1, 1; x) + 7H(0, 0, 0, 0; x) 

Z o D 

+ 10 H(l, 0, 0, 0; x) +4H(1, 1,0,0; x) + H(l, 1, 1, 0; x)) + O (e^)] , 

(A.14) 
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Ci234 =27riSr (Q') ' ^'{l-x) 



2 6 . 1 /57r2 



+ 34C3 + — vr^ H(0; x) 
3 



+ -TT^ H(l; x) + 8 H(0, 1, 0; x) - 36 H(0, 0, 0; x) + 12 H(l, 0, 0; x) + 2 H(l, 1, 0; x) 
o 

+ e (-^TT^ + 80 H(0; x) Ca - 22 H(l; x)C3 + ^tt" H(0, 1; x) + y tt^ H(0, 0; x) 

- ^7r2H(l,0;x) + ^7r2H(l,l;x) + 16H(0,0, l,0;x) - 8H(1, 0, 1, 0; x) 
+ 40 H(0, 1, 0, 0; x) + 8 H(0, 1, 1, 0; x) - 76 H(0, 0, 0, 0; x) + 32 H(l, 0, 0, 0; x) 
+ 4H(l,l,0,0;x) +2H(l,l,l,0;x)) + O (e^)] , 

(A.15) 



-,2^-3-2e 

{1 - X) 

3 9 , . 1 /ISvr^ 



Vi2345 =2niSr (Q") ' (1 - x)-^-^'x^ 



-- - - H(0; x) + - f — ^ + 3 H(0, 1; x) - 22 H(0, 0; x) + 2 H(l, 0; x)) + 23(3 

^^7r2H(0;x) + l-. 
-3H(0,l,l;x) -48H(0,0,0;x) +8H(l,0,0;x) + O (e) 



+ ^TT^ H(0; x) + ^TT^ H(l; x) + 10 H(0, 0, 1; x) - 2 H(l, 0, 1; x) + 10 H(0, 1, 0; x) 



(A.16) 



B. 0-DIS coefficient functions 



The coefficient functions for (/)-DIS structure functions derived by comparing them with 
the antennae as done in section 7 are given by: 



L Z Z o 



- z(4CaH(0; z) + 4CaH(1; z)) + z^{ACa^[^\ z) + 4CaH(1; z)) 



(B.l) 



-2CA + CA7r +(5(l-z)( ^ + + T^) + T 



27 



24 



12 



18 



^f^)do + (-^^ + + ^^"^2 - 4CaH(0, O; z) - 4CaH(0, 1; z) 

4CaH(1, O; ^) - 4CaH(1, 1; z) + (^CaH(0; z) - JiV^ H(0; z) + 2C^H(0, 0; z) 



+ 2C^H(0, 1; z) + 2CaH(1, 0; z)) + - ( 



1 6 " ' ' ' 3 
1,67Ca SiVp Catt^ 
2 ^ 18 9 



2 +^CaH(0;z) 
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- ^ArpH(0; z) + y CaH(1; z) - ^AT^ H(l; z) + 2CaH(0, 0; z) + 2CaH(0, 1; z) 

+ 2C^H(1, 0; z) + 2C^H(1, 1; z)) + z( 1— + 2CaH(0, 0; z) + 2C^H(0, 1; z) 

+ 2CaH(1, 0; z) + 2CaH(1, 1;z))+ - 2CaH(0, 0; z) - 2CaH(0, 1; z) 

-2CaH(1,0;z) -2CaH(1,1;2))1 , (B.2) 



4CfH(0; z) + 4CfH(1; z) - ^(3Cf + 4CfH(0; z) + 4CfH(1; z)) 

+ z(2Cir - 2CfH(0; z) - 2CirH(l; z))] , (B.3) 



-2Cf + — ^ 2CfH(0, 0; z) - 2(7fH(0, 1; z) - 2CfH(1, 0; z) - 2CfH(1, 1; z) 

+ \C-^- ^ + ^CfH(0; z) + ^CfH(1; z) + 2Ci.H(0, 0; z) + 2Cf H(0, 1; z) 

+ 2Cj.H(l, 0; z) + 2CfH(1, 1; z)) + z{ ^ CfH(0; z) - CirH(l; z) 

+ CfH(0, 0; z) + Cj.H(0, 1; z) + CfH(1, 0; z) + CfH(1, 1; z))] , (B.4) 



5150CI _ 566C.iV. 521C.iV. _ 266C|^ 20 _ 
27 27 9 9 9 

^5^^,28 2 4iV|7r^ 101C>^ 242 .. 28 

+ 68ClC3)do + (^Cl - ^-^CaN, + ^iV| - ^^CW)d, - (f - ^CAiV^)d2 

S2 449 242 

+ 8Cld3 + -CWHi-l; z) + — ClH(0; z) - — CAiVi.H(0; z) 

+ ^CFiVFH(0; z) - ^Clvr^HCO; z) - ^CfA^f^'H(0; z) + ii^ClH(l; z) 

o o o y 

256 62 40 532 

- — <7AiVpH(l; z) + -CfNfYL{1- z) - y CVH(1; z) + — C1h(0, 1; z) 

- ^CaA^fH(0, 1; z) + 12CfA^fH(0, 1; z) + 64C1h(0, O, 1; z) + 16CfA^fH(0, O, 1; z) 

-96C1H(0,-1,0;z) + 16CiH(-l, 0; z) - 32CiH(-l, 0, 1; z) + 202C1H(0, 0; z) 

- 12CyiAfFH(0, 0; z) + 14CFiVFH(0, 0; z) + 136C^H(1, 0; z) - 16CAiVFH(l, 0; z) 
+ 4CFiVFH(l, 0; z) + 136C1h(1, 1; z) - 16CAiVFH(l, 1; z) + 4CfA^fH(1, 1; z) 

+ 112C1h(1, 0, 1; z) + 64CiH(0, 1, 0; z) + 8CFAfFH(0, 1, 0; z) + 80C1h(0, 1, 1; z) 



+ 8CfNfR{0, 1, 1; z) + 64C1h(-1, -1,0; z) - 80C1h(-1, 0, 0; z) 

+ 20CfA^fH(0, 0, 0; z) + 80C\B{1, 0, 0; z) + 112C1H(1, 1, 0; z) + 96C1H(1, 1, 1; z) 

+ ^ (28C1C3 - y CVH(-1; + ^CVH(0; z) - 8C1h(0, 0, 1; z) 

+ 24C1h(0, -1, 0; z) + 16C1h(-1, 0, 1; z) - 32C1H(-1, -1, 0; z) 

1 778 

+ 40ClH(-l, 0, 0; z) - 12C1h(0, 0, 0; z)) + (- c1h(0; z) 

1 — z ,J 

+ ^CANFm z) - ^A^|H(0; z) + y Cl7r2H(0; z) - ^C1H(0, 1; z) 

+ —CaNfR{0, 1; z) - 56C1H(0, 0, 1; z) + 24C1H(0, -1, 0; z) - 44C1H(0, 0; z) 
3 

1 76 32 
+ SCANpHiO, 0; z) - — ClH(l, 0; z) + —CaNfII{1, 0; - 56C1h(1, 0, 1; z) 

- 48C1H(0, 1, 0; z) - 56C1H(0, 1, 1; z) - 28C1H(0, 0, 0; z) - 40CiH(l, 0, 0; z) 

- 50CiH,M,0;.)) + + + 551^ 

1 'i^fi 88 fi4 8 

- ^ClH(0; z) + —CANFil{0-, z) + -^CFiVFH(0; - -iV|H(0; z) + SCIt^^11{Q- z) 

S40 88 fi4 8 20 

- — ClH(l; z) + y C^iV^H(l; z) + -Ci.iV^H(l; - -iV2H(l; z) + y CVh(1; z) 

32 16 

- 88C1h(0, 1; z) + — CAiVFH(0, 1; z) + — CFiVFH(0, 1; z) - 48C1h(0, 0, 1; z) 

o o 
88 220 
+ —C\YL{-1, 0; z) - 16ClH(-l, 0, 1; z) - — C1H(0, 0; z) + SCaNfR{Q, 0; ^) 

+ —CfNfR{Q, O; z) - — C1h(1, 0; z) + —CaNfR{1, O; z) + —CfNfR{1, 0; 2) 

352 32 16 

- — ClH(l, 1; z) + y CAiVi.H(l, 1; z) + y CfAT^^HCI, 1; ^) - 56C1h(1, O, 1; ^) 

- 48C1h(0, 1, O; z) - 56C1h(0, 1, 1; z) + 32C1h(-1, -1, O; z) - 40C1h(-1, O, 0; z) 

- 16C1H(0, 0, 0; z) - 40ClH(l, 0, 0; z) - 56C1H(1, 1, 0; z) - 48CiH(l, 1, 1; z)) 

, 4486C| 580CUiVF 386CWVf ^ ISeC^Tr^ 16 2 ^ 2 

+ ^( 2^ + ^7 9 + ^ yC^iV^TT --CfNf-k 

+ 72ClC3 - SCFiV^Ca + y CItt^HI-I; z) - ^CiH(0; z) + i^C^iV^ H(0; z) 

- hjFNF^: z) + y CVhCO; z) - ^C^iV^.7r2H(0; z) - ^C1H(1; z) 

182 38 20 224 

+ — CAiVFH(l; 2) - y CFiV^H(l; z) + y CVh(1; - — C1h(0, 1; z) 

32 

+ — CaA^fH(0, 1; z) + 8CfA^fH(0, 1; z) - 112C1h(0, O, 1; z) + 16CfA^fH(0, O, 1; z) 

+ 16(7lH(-l, O; z) - 16ClH(-l, 0, 1; z) - 66ClH(0, 0; z) + 12CaA^fH(0, O; z) 
+ 10CirArpH(0, 0; z) - 92C^H(1, 0; z) + 8Ca-/VfH(1, O; z) - 4CfA^fH(1, 0; z) 

- 92C1h(1, 1; z) + 8CaA^fH(1, 1; ^) - 4CfA^fH(1, 1; z) - 56C1H(1, 0, 1; z) 
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- 80C1h(0, 1, 0; z) + SCpMpiliO, 1,0; z) - 88CiH(0, 1, 1; z) + 8CfNfR{0, 1, 1; z) 
+ 32C1H(-1, -1, 0; z) - 40C|H(-L 0, 0; z) - 96CiH(0, 0, 0; z) 

+ 20CfA^fH(0, 0, 0; z) - 40C1h(1, 0, 0; z) - 56C1h(1, 1, 0; z) - 48C1h(1, 1, 1; z)) 
, 2/5747C| 554CaNf mCFNp 2420^^ , V 2 ,8 o 

- 96C1C3 + y Cl7r2H(-l; z) + ^C1H(0; - il^C^iVf H(0; z) - H^C^ iV^H(0; 

- f CVH(0; z) + ^ClH(l; z) - ^C^AT^ H(l; z) - l^C^iV^H(l; z) 

- ^cVh(1; z) + ^ClH(0, 1; z) - SCANpHiO, 1; z) - y C^iVFH(0, 1; z) 

+ 64(7iH(0, 0, 1; z) - 48(7iH(0, -1, 0; z) + ^C1H(-1, 0; 2) - l6ClR{-l, 0, 1; z) 
+ — ClH(0, 0; z) - 8CAiVFH(0, 0; z) - —CpNEiliO, 0; z) + — C1h(1, 0; z) 

- SCANpHih 0; ^) - y CFiV^H(l, 0; z) + ^C1h(1, 1; z) - SC^iV^ H(l, 1; z) 

- ^CfNf11{1, 1; z) + 56C1h(1, 0, 1; z) + 48(7lH(0, 1, 0; 2) + 56C1h(0, 1, 1; z) 

+ 32C1h(-1,-1,0;z) - 40C1h(-1, 0, 0; z) + 40C1h(0, 0, 0; z) + 40C1h(1, 0, 0; z) 
+ 56C1h(1, 1,0; z)+ 48CiH(l, 1, 1; z)) 1 , (B.5) 



27 27 9 3 9 

8 8 20 

- UCaCfC,^ - 8CU3 + -CACFTT^Hi-l; z) + -CVh(-1; z) - y C^C^ H(0; z) 

+ 15C|^H(0;z) - ^CfNfR{0;z) + ^CACir7r2H(0; z) - ^CVh(0;z) 

242 llfi 20 

+ — CaC^H(1; z) + 56C2 H(l; z) - — CpiY^Hll; z) - —CACFn^m; z) 

+ ^CaCfH(0, 1; z) + 20C|.H(0, 1; z) - ^CfNf^{0, 1; z) - 16CaCfH(0, 0, 1; z) 

+ 32C|^H(0, 0, 1; z) - 16CaCfH(0, -1, 0; z) - 32C|^H(0, -1, 0; z) 
+ 32CACirH(-l, 0; z) - 16C|^H(-1, 0; z) - 16CaCfH(-1, 0, 1; z) 

+ — CaCfH(0, 0; z) - 12C|H(0, 0; z) - — CFiV^H(0, 0; z) + —CaCf^I, 0; z) 
+ 32C^H(1, 0; 2) - ^CfNf^{1, 0; z) + ^CaCf^{1, 1; ^) + 44C|H(1, 1; z) 

- \cfNf^{1, 1; z) + 40CaCfH(1, 0, 1; z) + 16C|H(1, 0, 1; z) + 24C|,H(0, 1, 0; z) 

+ 8CaCfH(0, 1, 1; z) + 24^2 H(0, 1, 1; z) + 32C|H(-1, -1, 0; z) 

- 24CaCfH(-1, 0, 0; z) - 16C|H(-1, 0, 0; z) - 40CaCfH(0, 0, 0; z) 

+ 20C|H(0, 0, 0; z) + 40CaCfH(1, 0, 0; z) + A{)CaCf^{1, 1, 0; z) + 16C|H(1, 1, 0; z) 
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+ 40C.C,H(1, 1. 1;.) + SCjHd, 1,1;.)+ '-(-'-^^^ - + 

+ I^CaCftt^ + eClvr^ - IcfNft^^ + 72CaCfC^ - 32C^C3 + aC Fi^^^i^i-l; z) 
y y o 

+ ^CVh(-1; z) - 1^CaCfH(0; z) - 10C|H(0; z) + ^CfNf^; z) 

1 fi 8 

+ — CaCf^^H(0; z) + -C|7r2H(0; z) - 32C^C7^H(1; z) - 64C7|H(1; z) 

o 

+ H^CpiVf H(l; z) + ^CaCf 7r2H(l; z) - 1^CaCj.H(0, 1; z) - SGC^ H(0, 1; z) 

1 fi 88 

+ — CfATfHCO, 1; z) - 32CaCfH(0, 0, 1; z) - 16C^H(0, 0, 1; z) + — CaCfH(-1, 0; z) 

o o 

212 32 

- 16CaCj.H(-1, 0, 1; z) - — CaCfH(0, 0; z) + y CfA^fH(0, 0; z) 

212 Ifi 

- —CaCf^{1, 0; z) - 36C|,H(1, 0; z) + — C^ATj^Htl, 0; z) - 56(7aCfH(1, 1; z) 

o o 



48C|.H(1, 1; z) + -CfNfB^{1, 1; 2) - 40CaCfH(1, 0, 1; z) - 16C|H(1, 0, 1 
o 



- 32CaCfH(0, 1, 0; z) - 16C|H(0, 1, 0; z) - 40CaCfH(0, 1, 1; z) - 16C|.H(0, 1, 1; z) 
+ 32C|H(-1,-1,0;z) - 24CaCfH(-1, 0, 0; z) - 16C|.H(-1, 0, 0; z) 

- 16CaCfH(0, 0, 0; z) - 40CaCfH(1, 0, 0; z) - 40CaCfH(1, 1, 0; z) 

- 16C|,H(1, 1, 0; z) - 40CaCfH(1, 1, 1; z) - SCf, H(l, 1, 1; z)) + z(^^^^ - 15(7^ 

+ - y CaCfvt^ + ^ - ^CfNf.^ + 44CaCK3 - 120^3 

+ \caCfi^'^^{-1; z) + ^C|7r2H(-l; z) - 57CaCfH(0; z) + 18C|,H(0; z) 

22 16 8 1 27 

+ —CfNf^{Q\ z) + -CaCft^M^; z) + -C^vr2H(0; z) - _CaCfH(1; z) 

+ 3C^H(1; z) + ^CfNf^; ^) + ^CACFTr'^m-. ^) + y CaCfH(0, 1; z) 

- 24C|H(0, 1; z) + -CfA^fH(0, 1; z) - 32CaCfH(0, 0, 1; z) - 16C|H(0, 0, 1; z) 

3 

26 

+ 12CaCfH(-1, 0; z) - 16C|.H(-1, 0; z) - 8CaCfH(-1, 0, 1; z) + — CaCfH(0, 0; z) 

+ C|.H(0, 0; z) + ^CFiVFH(0, 0; z) + ^CaCfH(1, 0; z) - 14C^H(1, 0; z) 

8 26 4 

+ -CFiVfHCl, 0; z) + y CaCfH(1, 1; z) - 24C^H(1, 1; z) + -CfA^fH(1, 1; z) 

- 20CaCfH(1, 0, 1; z) - 8C|^H(1, 0, 1; z) - 24C'aCfH(0, 1, 0; z) 

- 12C|.H(0, 1, 0; z) - 28CaCfH(0, 1, 1; z) - 12C|H(0, 1, 1; z) + 16C|^H(-1, -1, 0; z) 

- 12CaCfH(-1, 0, 0; z) - 8C|H(-1, 0, 0; z) - 28CaCfH(0, 0, 0; z) 

- 10C|.H(0, 0, 0; z) - 20CaCfH(1, 0, 0; z) - 20CaCfH(1, 1, 0; z) - 8C|.H(1, 1, 0; z) 

- 20CACFm 1, 1; - 4C|,H(1, 1, 1; z)) + ^2( 608CaCf _ ^CaCftt' 
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- — C^CiT H(0; z) - — C^CiT H(l; z) + 16CaCfH(0, 1; z) + —CaCf^^-^ 0; ^) 

64 16 16 1 

+ y CaCfH(0, 0; z) + -CACpm. 0; z) + -CACpm, 1; ^))J • (B.6) 

They are in full agreement with [56]. 
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